i     1L 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


T  HES I S 

<b>4  ii>Tr 


DIRECT  BIT  DETECTION  RECEIVER 

PERFORMANCE  ANALYSES 

FOR  8-DPSK  AND  16-DPSK  MODULATED  SIGNALS 

OPERATING 

WITH  IMPROPER  CARRIER  PHASE 

SYNCHRONIZATION 

by 

Mehmet  Sevki  Sekerefeli 


December  1987 
Thesis  Advisor        Daniel  C.  Bukofzer 


Approved  for  public  release;  distribution  is  unlimited. 


T239231 


'R'~v  C-ASS  :  CA~  ON  'J; 


S    =iG: 


REPORT  DOCUMENTATION  PAGE 


!EPOR*  SEC-R;~V  CLASSY  CATiON 

NCLASSlflED 


ID    RESTRICTIVE    MARKINGS 


ECuRlTY  CLASSIFICATION  AUThORiTY 


>ECLASSiFiCAT:ON  /  DOWNGRADING   SCHEDULE 


3.  DISTRIBUTION /AVAILABILITY  OF   REPORT 

Approved  for  public  release  ; 
distribution  is  unlimited. 


RFORMING  ORGANIZATION  REPORT  NuMBER(S) 


5    MONITORING  ORGANIZATION   REPORT   NUMBER(S) 


JAME  OF  PERFORMING  ORGANIZATION 

val    Postgraduate  School 


60    OFFICE  SYMBOL 
(if  applicable) 
62 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Naval    Postgraduate  School 


.ODRESS  (Gfy,  Srjfe,  and  ZiP  Code) 

nterey,  California  93943-5000 


7o.   ADDRESS  (City.  State,  and  ZIP  Code) 

Monterey,   California      93943-5000 


iAME  OF  FUNDING.  SPONSORING 
ORGANIZATION 


8d    OFFICE  SYMBOL 
(If  applicable) 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION   NUMBER 


.DDRESS  (Gfy,  State,  ana  ZIP  Code) 


:-0    SOURCE  Or  FUNDING   NUMBERS 


PROGRAM 
ELEMENT  NO. 


PROJECT 

NO. 


TASK 
NO 


WORK   UNIT 
ACCESSION  NO 


ITLE  {include  Security   Classification) 

DIRECT  BIT  DETECTION  RECEIVER   PERFORMANCE  ANALYSES   FOR  8-DPSK  AND   16-DPSK  MODULATED 
3NALS   OPERATING   WITH   IMPROPER  CARRIER    PHASE  SYNCHRONIZATION. 


ERSONAi  A'JTHOR(S) 

Sekerefeli,  Mehmet  Sevki 


TYPE  O-  RE°Oa7 

Master's  Thesis 


3b    TiME  COVERED 
=ROM  -O 


14    DATE  0=  RE°ORT    (Year,  Month,  Day) 

1987  December 


PAGE   COUNT 
131 


UPOLEVENTAR*  NOTATION 


:  ELD 


tOl  ° 


"6    So3>£CT  TERMS  {Continue  on  reverse  it  necessary  and  identity  by  block   number) 

Communication,   Digital    Communication, DPSK,   PSK,  M-ary  DPSK, 
Direct  Bit  Detection,   Modulation,   Demodulation. 


3STRACT  '.Continue  on  reverse  it  necessary  anrj  laentity   Dy  blocK  number) 

The  application  of  Direct  Bit  Detection  methods  are  analyzed 
laluated  in  digital  communication  systems  employing  Differential  F 
lift  Keyed  (DPSK)  modulaticn.  Assuming  an  additive  white  Gaus 
)ise  interference  model,  and  direct  bit  detection  methods  used  in 
Iceiver,  S-DPSK  and  16-DPSK  communication  systems  are  considered 
leir  performances  evaluated  in  terms  of  the  delivered  bit  error  rat 
function  of  signal-to-noise  ratio.  The  advantages  and  disadvant 
I  Direct  Bit  Detection  Receiver  (DBDR)  systems  used  in  conduction 
Inferential  phase  encoding  are  determined,  with  specific  applicatio 
•DPSK  and  16-DPSK  modulated  signals.  The  effect  of  improper  rece 
irrier  phase  synchronization  is  considered  in  detail,  and  resul 
irformance  degradations  are  evaluated.  numerical  results  show  that 
iceiver  phase  errors  of  more  than  a  few  degrees,  severe  oerform 


and 
hase 
sian 

the 

and 
e  as 
ages 
with 
n  to 
iver 
ting 

for 
ance 


nSTR'8uT:0\    AVAlLABlLll  v  OF   A3STRACT 
UNCLASSIF  E3  UNLIMITED       D   SAME   AS  RPT  fj  QTiC  USERS 


21     ABSTRACT  SECURITY   CLASSIFICATION 

imriASSTFTFn 


IAIVE  O-  sES-'0'.S  Bl£   !•«:>■«  D    A: 

Professor  D.   C.    Bukofzer 


22b  TELEPHONE  (Include  AreaCode) 

(408)    645-2859 


2c    OFFICE  SYMBOL 

62Bh 


ORM  1473.  34  mar 


33  ADR  edition  may  oe  useo  untn  exnausted 
All  other  editions  are  obsolete 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

O  U.S.  Gov«rnm«nt  Printing  C'lce.   1Si6 — 60S-24. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fTFh«n  Dmtm  Bnfwd) 


19  (Continued) 

degradations  result  for  both  DPSK  modulation  schemes,  unless  a 
complete  phase  reversal  (i.e. ,180°)  takes  place. 


S    N  0102-  LP-  014-  6601 


SECURITY   CLASSIFICATION  OF  THIS  F«AGE(T»»»»n  Dat»  Bnfrmd) 


Approved  for  public  release;  distribution  is  unlimited. 


Direct  Bit  Detection  Receiver  Performance  Analvses 

for  S-DPSK  and  16-DPSK  Modulated  Signals  Operating 

With  Improper  Carrier  Phase  Synchronization 


by 


Mehmet  Sevki  Sekerefeli 

Lieutenant  Junior  Grade,  Turkish  Navy 

B.S.,  Turkish  Naval  Academy,  1981 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  ELECTRICAL  ENGINEERING 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 

December  1987 


C( 


ABSTRACT 


The  application  of  Direct  Bit  Detection  methods  are  analyzed  and  evaluated  in 
digital  communication  systems  employing  Differential  Phase  Shift  Keyed  (DPSK) 
modulation.  Assuming  an  additive  white  Gaussian  noise  interference  model,  and  direct 
bit  detection  methods  used  in  the  receiver,  8-DPSK  and  16-DPSK  communication 
systems  are  considered  and  their  performances  evaluated  in  terms  of  the  delivered  bit 
error  rate  as  a  function  of  signal-to-noise  ratio.  The  advantages  and  disadvantages  of 
Direct  Bit  Detection  Receiver  (DBDR)  systems  used  in  conjuction  with  differential 
phase  encoding  are  determined,  with  specific  application  to  8-DPSK  and  16-DPSK 
modulated  signals.  The  effect  of  improper  receiver  carrier  phase  synchronization  is 
considered  in  detail,  and  resulting  performance  degradations  are  evaluated.  Numerical 
results  show  that  for  receiver  phase  errors  of  more  than  a  few  degrees,  severe 
performance  degradations  result  for  both  DPSK  modulation  schemes,  unless  a 
complete  phase  reversal  (i.e.,  180°)  takes  place. 
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I.  INTRODUCTION 

M-ary  differential  phase  shift  keyed  (M-DPSK)  modulation  is  a  widely  used 
technique  in  digital  communication  applications.  Bandwidth  efficiency,  relatively  good 
noise  immunity,  constant  signal  envelope,  and  simplicity  of  implementations,  make  this 
scheme  particularly  attractive  for  use  over  satellite,  terrestrial  radio  and  voiceband 
telephone  channels.  While  system  analyses  pertaining  to  the  performance  of  M-DPSK 
receivers  abound  in  the  literature,  treatment  is  usually  restricted  to  the  case  where 
signal  transmission  takes  place  over  an  additive  white  Gaussian  noise  channel.   [Ref.  1] 

Since  the  digital  information  transmission  is  based  on  carriers  modulated  by  the 
symbol  waveforms,  the  symbol  error  probabilities  which  quantify  receiver  performance 
can  often  be  computed  directly  as  demonstrated  by  the  many  derivations  of  such 
probabilities  found  in  the  literature  (  see  [Ref.  2:  pp.  228-234]  and  [Ref.  3:  pp.  200-212] 
for  example  ).  However,  from  the  point  of  view  of  the  recipient  of  the  digital 
information,  it  is  the  Bit  Error  Ratio  (BER)  that  becomes  important  in  digital 
communication  applications.  While  results  on  the  BER  performance  of  M-ary 
modulation  receivers  is  far  more  limited,  closed  form  expressions  for  the  BER  of  M- 
PSK  receivers  have  been  obtained  by  Lee  [Ref.  4],  with  similar  results  independently 
derived  by  Tan  [Ref.  5]. 

This  thesis  is  devoted  to  analyzing  the  application  of  direct  bit  detection  methods 
which  allow  the  recovery  of  bits  in  the  received  symbol  individually,  in  differential 
PSK  digital  communication  systems.  For  analysis  purposes,  8-DPSK  and  16-DPSK 
systems  are  considered  and  their  performances  evaluated.  The  particular  question  to 
be  answered  is  whether  for  8-DPSK  and  16-DPSK  modulated  signals,  direct  bit 
detection  methods  provide  comparable  performance  to  conventional  phase  detection 
receivers,  especially  under  conditions  of  improper  receiver  carrier  phase 
synchronization. 

In  this  introductory  chapter,  background  information  on  M-PSK  communication 
systems  is  given.  Known  optimum  receiver  structures  are  presented,  and  Direct  Bit 
Detection  Receivers  (DBDR)  are  introduced.  In  Chapters  II  and  III  a  presentation  of 
the  transmitter  and  receiver  structures  respectively  is  given  along  with  their  operational 
properties,   and  the  logic  implemented  in  the   DBDR  for   8-DPSK   and   16-DPSK 
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modulated  signals  is  discussed  in  detail.  Chapter  IV  includes  the  noise  performance 
analyses  of  8-PSK  and  8-DPSK  receivers,  and  Chapter  V  presents  similar  receiver 
analyses  for  16-PSK  and  16-DPSK.  In  both  of  these  chapters,  receiver  local  oscillator 
phase  errors  have  been  incorporated  so  as  to  be  able  to  account  for  degrading  effects 
due  to  improper  carrier  phase  synchronization.  The  obtained  mathematical  expressions 
of  receiver  performance  are  evaluated  on  a  computer  and  graphically  presented  as 
plots  of  bit  error  ratio  as  a  function  of  signal  to  noise  ratio  (SNR)  and  phase  errors. 

A.  BACKGROUND 

Direct  bit  detection  methods  utilize  two  signal  processing  channels  in  the  receiver 
in  order  to  recover  the  digital  information  without  phase  angle  measurements. 
Moreover,  these  methods  exhibit  optimum  performance  in  a  BER  sense  when  carrier 
synchronization  is  completely  achieved.  Because  of  the  attractiveness  and  simplicity  of 
DBDRs,  some  studies  have  been  carried  out  in  this  area,  and  a  modified  two-channel 
receiver  for  8-PSK  has  been  built  and  its  noise  performance  measured  by  Thompson 
[Ref.  6].  Furthermore,  noise  performance  evaluations  of  DBDR  for  8-PSK  have  been 
presented  by  Myers  and  Bukofzer  [Ref.  7],  while  similar  analyses  for  the  16-PSK  case 
can  be  found  in  Reference  8.  Advantages  of  direct  bit  detection  methods  are  well 
documented  in  these  referenced  works,  in  which  all  analyses  have  been  done  assuming 
no  receiver  local  oscillator  phase  error,  Gray  Code  bit  to  symbol  mapping  and  an 
additive  white  Gaussian  noise  interference  model.  We  present  here  8-PSK,  8-DPSK, 
16-PSK  and  16-DPSK  communication  systems  under  similar  assumptions,  except  now 
the  receiver  has  an  assumed  known  local  oscillator  phase  error. 

B.  THEORY  OF  M-PSK  COMMUNICATION  SYSTEMS 
1.  Representation  of  M-PSK  Signals 

For     M-PSK     modulation     with     equal     signal     energies,     a     convenient 
representation  [Refs.  3:  p.p.  192-205]  of  the  signal  set  is  given  by 

s.(t)  =    V2E~Tscos[27tfot  +  27t(i-l)/M],    i=  1,2,...,M     ,    0  <  t  <  Tg  (1.1) 

T 
s 

Es=  J    S;2(t)dt   ,  i=l,2,..,8  (1.2) 

0 

Ts  =  symbol  duration 

fQ  =  carrier  frequency  (Hz.) 
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for  convenience  f  >  >   1/T  is  assumed, 
o  s 

A  suitable  orthonormal  signal  set  for  the  representation  of  the  above  signals  is 


^(0=  V2/Ts  cos(27tfot) 

q>2(t)  =  v'2/Ts  sin(27tfot)         ,      0  <  t  <  Ts  (1.3) 

so  that  by  using  trigonometric  identities,  s}  (t)  can  be  expanded  in  terms  of  cpj(t)  and 
<p,(t).   That  is, 

s.(t)  =    (p^t)  JWS  cos[  27C(i-l)/M  ]  -  <p  2(t)  V2E;  sin[  2w(i-l)/M  ]  (1.4) 

where    i=  1,2 M  . 

A  plot  of  the  signal  constellation  for  M-PSK  is  shown  in  Figure  1.1  for 
various  values  of  M,  where  the  coordinates  of  the  i  signal  vector  representing  the 
signal  s.(t)  are 

{V2ES  cos[2ir(i-l)/M] }   ,-   {  V2ES   sin[27t(i-l)/M] }  (1.5) 

For  M-DPSK,  the  signal  representations  are  exactly  same  as  those  already  presented, 
except  that  the  digital  information  is  transmitted  in  terms  of  phase  differences  between 
consecutive  signal  transmissions,  rather  than  absolute  signal  phases. 
2.  Gray  Coding  for  M-PSK  Signals 

Before  proceeding  with  the  Gray  coding  for  M-PSK  signals,  we  introduce 
some  definitions  which  will  henceforth  be  used 

1.  Symbol  :  Transmitted  or  received  signal  in  sinusoidal  form  (s}(t)  where 
i=  1,2,...M)  which  is  related  to  the  symbol  state. 

2.  Symbol  state  :  The  binary  digit  assignment  to  the  symbols  which  always 
consists  of  a  block  of  log2M  bits  in  length^  where  i=  1,2,.. .,M.) 

3.  Data  state  :  Block  of  log2M  bits  corresponding  to  the  data  that  must  be 
delivered  to  the  intended  user,  and  which  is  processed  at  the  encoder  input  to 
determine  the  symbol  to  be  transmitted.   (D.  where  i=  1,2,. ...,M.) 

The  relationship  between  symbol  error  probability  and  bit  error  probability  in 

the  corresponding  log2M  bit  groups  depends  upon  the  assignment  of  bits  to  the  symbol 

states.    The  preferred  bit  to  symbol  state  mapping  in  most  cases  is  the  Gray  Code,  in 

which  the  symbol  states  which  correspond  to  adjacent  phases  differ  in  only  one  bit 

position.    This  is  because  in  the  coherent  demodulation  process  the  most  likely  error 
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Figure  1.1     Signal  Constellations  lor  (a)  BPSK  (b)  QPSK  (M  =  4)  (c)  8-PSK  (d)  16-PSK. 
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involves  mistaking  a  correct  symbol  with  its  nearest  neighbor  and  so  such  symbol 
errors  translate  into  single  bit  errors. 

In  Direct  Bit  Detection  Receivers  (DBDR),  Gray  code  assignment  to  the 
symbol  states  are  selected  for  the  above  given  reason.  In  this  thesis,  a  modified  Gray 
code  is  used,  which  for  comparison  purposes  is  shown    along  with  the  conventional 


TABLE  1 

COMPARISON  OF  CONVENTIONAL  AND  MODIFIED  GRAY 

CODING  SCHEMES 

Svmbol 
States 

Grav 

Coding 

Modified  Gray 
Coding 

8-PSK 

16-PSK 

8-PSK 

16-PSK 

s, 

000 

0000 

111 

mi 

s2 

001 

0001 

101 

1101 

2} 

011 

0011 

001 

1001 

s4 

010 

0010 

011 

1011 

s5 

110 

0110 

010 

0011 

56 

111 

0111 

000 

0001 

2, 

101 

0101 

100 

0101 

Ss 

100 

0100 

no 

0111 

s9 

1100 

0110 

Sio 

1101 

0100 

£11 

mi 

0000 

«U 

1110 

0010 

Sn 

1010 

1010 

«14 

1011 

1000 

fin 

1001 

1100 

s,« 

1000 

1110 

Gray  Code  in  Table  1. 

Using  the  modified  Gray  Code  on  Table  1,  the  vectors  representing  the 
symbols  with  components  along  the  q>1  and  <p2  axes  have  been  shown  (see  Figure  1.2) 
such  that  r,    is  the   symbol  component    along  the  <Pj    axis  and  r2  is  the  symbol 
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component  along  the  <p2  axis.  Therefore  for  M  =  8,  (see  Figure  1.2.a)  in  order  to 
recover  Most  Significant  Bit  (LB)  of  the  transmitted  symbols  correctly,  rj  >  0  whenever 
the  symbol  states  g,,  g,,  S?  and  S8  are  transmitted  and  r}  <0  whenever  the  symbol 
states  g3,  54,  55  and  S6  are  transmitted.  The  Middle  Bit  (MB)  of  the  transmitted 
symbols  can  be  recovered  correctly  if  (r^2  -  (r2)2  >  0  whenever  the  symbol  states  S{, 
54,  S5  and  5g  are  transmitted  and  (v^)2  -  (r2)2  <  0  whenever  the  symbol  states  S2,  £3, 
S6  and  S7  are  transmitted.  The  Least  Significant  Bit  (RB)  is  recovered  correctly  if 
r2>0  whenever  the  symbol  states  g,,  S?,  S3,  and  g4  are  transmitted  and  r2<0 
whenever  the  symbol  states  55,  S6,  g7,  and  Sg  are  transmitted. 

Similar  logic  has  been  set  up  for  the  M=  16  case  (see  Figure  1.2.b)  but  now 
each  symbol  is  related  to  4  bits  in  a  symbol  state.  Therefore  for  this  case,  4  bits  must 
be  recovered  from  each  received  ^7TV,bol,  where  the  bits  in  a  symbol  state  are  labeled 
Most  Significant  Bit  (LB),  Middle  Bit  Right  (MBR),  Middle  Bit  Left  (MBL),  and  Least 
Significant  Bit  (RB).  For  example  from  Figure  1.2.b,  S7  is  0101,  corresponding  to 
LB  =  0,  MBL=1,  MBR  =  0  and  RB=1.  In  order  to  recover  the  LB  correctly,  it  is 
necessary  that  ^>0  for  transmitted  symbol  states  g,,  S2,  g3,  54,  513,  S14,  £>15,  and 
516  and  that  v^<  0  for  transmitted  symbol  states  S5,  56,  S7,  S8,  S9,  S1Q,  g,,,  and  SJ2. 
If  r2>0  whenever  the  symbol  states  g,,  52,  S3,  S4,  S5,  56,  S7,  and  58  are  transmitted 
and  r2<0  whenever  the  symbol  states  59,  S10,  g,,,  512,  5]3,  Si4,  S15,  and  S16  are 
transmitted,  the  RB  will  be  recovered  correctly.  If  (r^2>  (r2)2  for  transmission  of 
symbol  states  5p  S2,  57,  S8,  5g,  510,  515,  and  S16,  and  (rj)2<(r2)2  for  transmission  of 
symbol  states  53,  54,  55,  S6,  5n,  S12,  S13,  and  514,  the  MBL  will  be  recovered 
correctly.  Finally,  if  [  (r^2  -  (r2)2  ]2  -  (2  r2r2)2  >  0  upon  transmission  of  the  symbol 
states  Sp  54  ,  S$,  S8,  Sg,  Sn,  S{V  and  S{6,  and  [  (rL)2  -  (r2)2  ]2  -  (2  r^)2  <  0  upon 
transmission  of  the  symbol  states  on  the  states  S2,  53,  S6,  S7,  S10,  g,,,  S14,  and  g15, 
then  the  MBR  will  be  recovered  correctly. 

In  subsequent  chapters,  the  receiver  logic  and  its  performance  will  be 
explained  in  detail.  Additionally,  application  of  this  methodology  will  be  presented  in 
so  far  as  differential  PSK  modulation  is  concerned.  Before  doing  so  however,  a 
summary'  of  known  results  on  the  performance  of  conventional  and  DBDR's  for  M- 
PSK  and  M-DPSK  is  presented  next. 

3.  Known  Results  on  the  Performance  of  M-PSK  and  M-DPSK  Receivers 

In  this  section,  known  results  on  the  receiver  performance  for  M-PSK  and 
DPSK  modulated  signals  are  presented  for  two  cases,  namely  conventional  receivers 
and  Direct  Bit  Detection  Receivers  (DBDR). 
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Figure  1.2    Symbol  State  Assignments  for  (a)  8-PSK  (b)  16-PSK. 
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a.  Conventional  M-PSK  and  M-DPSK  Receivers 

Coherent  demodulation  of  M-PSK  implies  in  principle  the  use  of  M  signal 
processing  channels  in  the  receiver,  as  shown  in  Figure  1.3  (from  Ha  [Ref.  10]  )  where 
different  forms  of  the  optimum  receiver  structure  of  M-PSK  signaling  have  been  given. 
Performance  analysis  results  for  optimum  receivers  are  well-documented  in  terms  o[ 
symbol  error  rate  as  a  function  of  signal  to  noise  ratio  [Refs.  2,9],  while  more  recent 
analyses  have  focused  on  the  receiver  probability  of  bit  error  (see  [Ref.  4]  and 
[Ref.  3:  p.p.  204-207  ]  for  example).  Symbol  error  probabilities  can  be  computed  from 
knowledge  of  the  channel  and  signal  characteristics.  However  for  transmission  of 
binary  data,  when  analyzing  system  performances  as  a  function  of  different  levels  of 
modulation,  the  bit  error  probability  rather  than  the  symbol  error  probability  is  of 
interest,  as  previously  explained.  In  Prabhu  [Ref.  1:  p.  198],  the  bit  error  prrbJ>ility 
has  been  presented  only  for  the  case  of  M-ary  orthogonal  signal  sets.  However,  since 
the  M-PSK  signal  set  is  not  orthogonal,  those  results  are  not  applicable  here.  The 
closed  form  expression  for  the  BER  of  M-PSK  has  been  determined,  when  a  Gray 
Code  bit  mapping  is  used  [Ref.2:  p.  198].   The  results  can  be  summarized  as  follows 

n 
Pb(M)  =  £  (k/n)  Pk  (1.6) 

k=l 

where 

M  =  Number  of  distinct  signal  waveform, 

Pk  =  Probability  of  k  bits  in  error  in  a  received  n  bit  data 

block  [Ref.  5:  pp.24-32], 
n  =  log2  M  , 
Pb(M)  =  Bit  error  probability. 

A  derivation  of  Equation  1.6  and  its  numerical  evaluation  (see  Figure  1.4  )  has  been 
carried  out  by  Tan  [Ref.  5:  p.33].  For  M-PSK  signalling,  receiver  performance  analyses 
in  terms  of  bit  error  rate  have  also  been  presented  by  Lee  [Ref.  4:  p.491],  and  plotted  in 
Figure  1.5. 

b.  Direct  Bit  Detection  Receivers  (DBDR's) 

DBDR  structures  for  M-PSK  have  been  derived  for  various  values  of  M 
[Refs.  6,7,8].  As  previously  indicated,  this  method  is  particularly  well  suited  for  M-PSK 
communication  systems.   The  DBDR  structures  for  8-PSK  (see  Figure  1.6)  and  16-PSK 
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Figure  1.3     Block  Diagram  of  Optimum  Receiver  (a)  Correlator  Realization 
(b)  Matched  Filter  Realization  (from  Ha  [Rcf.  10]). 
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Figure  1.4     Bit  Error  Probabilities  Versus  E^Nq  for  16-PSK,  8-PSK  and  QPSK  [Ref.  5  p.3l]. 


20 


n 


^5- 

PQ 


-5 


\ 


LEGEND 
□  8-PSK 
o     16-PSK 


T 


^ 


5  10 

SNR  (DB) 


-ty 


15 


20 


Figure  1.5     Probability  of  Bit  error  Versus  Eb/NQ  [Ref.4]. 
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(see  Figure  1.8)  have  been  presented  and    their  performances  have  been  derived  in 
References  6,7,8.   The  BER  is  quantified  by  the  probabilities  of  a  bit  error,  denoted  PE. 
For  8-PSK,  PE  has  been  derived  by  Myers  and  Bukofzer  [Ref.  7:  p.  13] 
where 


PE=  1/3  {  Q[  2  VSNR  cosa  ]  +    Q[  2  VSNR  sina  ]  + 


+  Q[  2  VSNR  (  cosa  +  sina)  ]    +  Q[  2  VSNR  (  cosa-  sina)  ]  - 


-  2  Q[  2  VSNR  (  cosa  -  sina)  ]  .  Q[  2  VSNR  (  cosa+  sina)  ]  }  (1.7) 

where 

SNR=Es/2N0=3Eb/2N0 

Es  =  Symbol  energy 

N«=  Noise  power  spectral  density  level  (one  sided) 

a    =  Basic  signal  phase  angle  (see  Figure  1.2). 


and 


CO 

,2 


Q(x)  =  J    [  1/V2TC  ]   exp(-u2/2)  du  (1.8) 


The  derivation  of  Equation  1.7  assumes  an  AWGN  channel,  Gray  coded  symbols  and 
equal  a  priori  bit  probabilities  (i.e.  ,  1/2).  Receiver  performance  in  BER  sense  is 
plotted  in  Figure  1.7  for  various  values  of  a  as  a  function  of  SNR. 

For  the  16-PSK  case,  probability  of  bit  error  is  derived  by  Bukofzer 
[Ref.  8]  and  the  resulting  receiver  performance  is  plotted  in  Figure  1.9  for  various 
values  of  the  angles  a  and  P,  where  P  is  shown  on  Figure  1.2.b. 

In  both  analyses,  it  has  been  assumed  that  the  receiver's  local  oscillator 
does  not  have  any  phase  error.  In  this  research,  similar  analyses  to  those  conducted 
in  References  7  and  8  have  been  carried  out  by  considering  the  effect  of  a  receiver 
local  oscillator  phase  error  for  8-PSK  and  16-PSK  cases,  as  well  as  the  noise 
performance  of  8-DPSK  and  16-DPSK  DBDR,  presented  as  BER  as  a  function  of 
SNR. 
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Fiaure  1.6    Block  Diasram  for  8-PSK  DBDR. 
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Figure  1.7    PE  Versus  SNR  (dB)  for  8-PSK  for  Various  Values  of  a. 
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Figure  1.8     Block  Diagram  for  16-PSK  DBDR. 
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Figure  1.9     PE  Versus  SNR  (dB)  for  Various  a  Values  for  16-PSK  DBDR. 
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II.  TRANSMITTER 

The  generation  of  differential  PSK  signals  can  be  performed  by  pre-processing 
the  data  bits.  This  chapter  is  devoted  to  analyzing  the  M-DPSK  transmitter  structure. 
One  general  form  of  the  transmitter  in  block  diagram  form  for  DPSK  communication 
systems  is  shown  in  Figure  2.1.  Explanation  of  the  general  elements  in  Figure  2.1  will 
be  included  along  with  the  derivation  of  the  transmitter  structure  for  8-DPSK  and 
16-DPSK  modulated  sienals. 


Output 

Input 

Conuorter 

Encoder 

Cam  ier 
Modu 1  at  1  on 

to 
Channa  1 

Do  lay 

Fisure  2.1     DPSK  Transmitter  Block  Diagram. 


A.       TRANSMITTER  STRUCTURE  FOR  8-DPSK  AND  16-DPSK  MODULATED 
SIGNALS 

The  main  task  in  generating  the  DPSK  transmitter  involves  the  design  of  the 

encoder,  which  depends  on  the  algorithm  used  to  relate  data  states  to  the  symbol  states 

prior  to  transmission  as  a  sinusoidal  waveform.    Basically,  the  transmitted  symbol  is 

the  quantized  equivalent  of  the  angle  between  the  successively  transmitted  symbols. 

For  reasons  previously  explained,  a  Gray  coding  scheme  is  utilized  in  the  symbol  state 

to  symbol  mapping  for  the  two  special  cases  of  M-DPSK  signalling    considered  here, 

namelv  M  =  8  and  M  =  1 6. 
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The  symbol  differential  phase  relative  to  the  data  state  to  be  transmitted  is  given 
on  Table  2  for  8-DPSK  and  on  Table  3  for  16-DPSK  modulated  signals.  In  both 
Tables,  it  is  assumed  that  an  equiphase  arrangement  for  the  transmitted  symbols  is 
used. 


TABLE  2 

GRAY  CODE  REPRESENTATION  OF  PHASES  FOR  8-DPSK  SIGNAL 

SET 

Differential 

LB 

MB 

RB 

Phase 

0 

0 

0 

0 

0 

0 

1 

a 

0 

1 

1 

2a 

0 

1 

0 

3a 

1 

1 

0 

4a 

1 

1 

1 

5a 

1 

0 

1 

6a 

1 

0 

0 

7a 

As  mentioned  earlier,  the  main  purpose  of  the  encoder  is  to  generate  a  mapping 
algorithm  whereby  the  phase  difference  of  consecutive  symbol  outputs  is  related  to  the 
input  data  in  accordance  with  the  information  in  Table  2  and  Table  3.  In  order  to 
achieve  successive  signal  transmissions  the  following  algorithm  is  implemented  in  the 
transmitter 

1.  Store  the  present  data  state  {D.,  i=  1,2...,M)  and  the  previous  symbol  state 
(Sj,i=l,2,3..,M). 

2.  From  the  Table  2  or  Table  3,  find  the  corresponding  differential  phase  value 
associated  with  the  present  data  state  as  Na,  where  N  is  an  integer  with  values 
ofO<N<M-l. 

3.  Generate  a  symbol  state,  such  that  the  corresponding  symbol  will  have  a  phase 
angle  Na  degrees  larger  (  in  additive  sense,  modulo  2tt)  than  the  previously 
transmitted  symbol. 
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TABLE  3 

GRAY  CODE  REPRESENTATION  OF  PHASES  FOR  16-DPSK  SIGNAL 

SET 

LB 

MBL 

MBR 

RB 

Differential 
Phase 

0 

0 

0 

0 

0 

0 

0 

0 

1 

a 

0 

0 

1 

1 

2a 

0 

0 

1 

0 

3a 

0 

1 

0 

4a 

0 

1 

1 

5a 

0 

0 

1 

6a 

0 

0 

0 

7a 

0 

0 

8a 

0 

1 

9a 

1 

1 

10a 

1 

0 

11a 

0 

1 

0 

12a 

0 

1 

1 

13a 

0 

0 

1 

14a 

0 

0 

0 

15a 

For  example,  assuming  the  previous  symbol  state  is  S7  (0101)  and  present  data  state  is 
D9  (0110)  which  corresponds  to  the  differential  phase  value  of  4a  (see  Table  3.).  From 
Figure  1.2.b,  proceed  4a  degrres  counterclockwise  from  g7,  resulting  in  the  symbol 
state  is  S14  (1111)  will  be  the  symbol  output  from  the  encoder. 

This  algorithm  can  be  shown  in  matrix  structure,  in  which  the  row  entries  and 
the  column  entries  of  the  matrix  are  the  symbol  states  and  data  states  respectively,  and 
the  elements  of  the  matrix  are  the  symbol  state  which  correspond  to  the  appropriate 
phase  difference  in  the  transmitted  symbol  (see  Figure   2.3   and  Figure  2.2).     For 
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example,  under  the  same  assumptions  given  in  the  previous  paragraph,  by  using  Q9 
(0110)  as  a  row  entry  and  5-  as  a  column  entry  in  the  matrix  in  Figure  2.3,  present 
symbol  output  from  encoder  is  found  as  (1111)  at  the  intersection  of  the  row  and 
column  of  the  matrix. 
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Figure  2.2    Encoder  Matrix  for  8-DPSK  Signals. 

The  hardware  implementation  of  the  encoder  structure  can  be  accomplished  in 
several  forms.  A  practical  approach  to  building  the  encoder  can  be  achieved  by  the  use 
of  EPROMs  (Erasable  Programable  Read  Only  Memory).  Once  the  truthe  table  the 
programming  of  the  EPROM  is  a  trivial  task  to  accomplish.  Demultiplexers  can  also 
be  used  to  build  the  encoder  circuitry;  however,  because  such  a  device  would  be  hard 
wired,  it  would  not  have  the  reprogrammability  of  the  EPROM  implementation. 

In  this  research,  the  8-DPSK  and  16-PSK  communication  systems  are  evaluated 
using  a  computer.  Appendix  A  includes  the  source  file  which  is  used  to  generate  the 
possible  outputs  of  an  M-PSK  and  an  M-DPSK  communication  system. 

The  carrier  modulator  uses  the  encoder  outputs  (symbol  states)  as  an  input,  and 
generates  the  symbols  (signals)  suitable  for  transmission.  The  signal  transmission 
involves  sinusoidal  voltages  whose  phase  is  related  to   the  symbol  states  (encoder 
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Fisure  2.3     Encoder  Matrix  for  16-DPSK  Signals. 
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outputs)  and  of  carrier  frequency  suitable  for  the  transmission  channel.  The  signal 
generation  can  be  achieved  via  a  multiplexer,  whose  structure  is  shown  in  Figure  2.4. 
For  this  kind  of  implementation,  (M-l)  phase  shift  elements  are  needed  with  set  phase 
shift  values.  Another  approach  would  involve  an  implementation  using  a  D/A 
converter  along  with  a  phase  modulator  which  can  be  obtained  with  a  VCO  (see 
[Ref.  11:  pp.269-271]. 


V2E,/Ts  cos[27tf0t] 


l  at.    bit 


aynnbo  I     <^ 
state 


2  nd.    bit 
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Figure  2.4    Carrier  Modulator  Structure  for  M-DPSK  Sisnalinu. 


B.       COMPUTER  GENERATION  OF  THE  TRANSMITrER  OUTPUTS 

For  computer  evaluation  of  the  M-ary  DPSK  communication  system  ,  a 
FORTRAN  program  called  BELIZ  has  been  written  (see  Appendix  A).  The  program 
generates  all  possible  data  state  and  symbol  state  combinations  for  M-DPSK,  and 
produces  the  corresponding  symbol  states  outputs.  The  program  outputs  are  given  in 
Appendix  B  in  tabular  form  along  with  the  transmitter  matrix  of  Figure  2.3  and  Figure 
2.2,  specified  in  decimal  form. 
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III.  RECEIVER 

A.       RECEIVER    STRUCTURE    FOR    8-DPSK    AND     16-DPSK    MODULATED 
SIGNALS 

The  introductory  chapter  presents  an  overview  of  the  structure  and  performance 

of  optimum  receivers  for  M-ary  signalling.    The  focus  of  this  chapter  is  on  DBDR's. 

As  mentioned  in  Chapter  II,  the  symbol  state  generation  based  on  the  data  state  and 

the  previous  symbol  state  is  the  most  important  task  performed  by  the  encoder. 

Therefore  in  the  receiver's  decoder,  a  similar  algorithm  must  be  performed  in  reversed 

sense,  that  is,  from  two  successive  symbol  states,  the  corresponding  data  state  must  be 

estimated.   A  generic  form  of  the  receiver  in  block  diagram  structure  is  shown  in  Figure 

3.1. 


User 

Rece i ued 

Demcdu lator 
and 

Carrier 

Recovery 

Decoder 

D/fl 

Converter 

data 

' 

1 

Ee  lay 

Figure  3.1     General  receiver  block  diagram  for  DPSK  signals. 

The  algorithm  for  recovering  the  transmitted  data  is  as  follows 

1.  Recover  the  symbol  state  from  the  received  symbol  using  direct  bit  detection 
techniques. 

2.  Store  the  past  symbol  state  5-,i=  1,2..., M. 

3.  Compare  the  present  symbol  state  with  the  past  symbol  state  to  determine  the 
incremental  "  phase  angle  "  between  the  two. 

4.  From  Table  2   or   3     find    the  corresponding  data  state  that  gave  rise 

to   the   "phase  angle"  between  the  two  consecutive  symbol  states. 
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For  example,  assuming  the  past  symbol  state  is  (56)  0001  and  the  present  symbol  state 
is  S3  (1001)  from  Figure  1.2.b,  the  angle  between  S3  and  S6  is  3a.  From  Table  3.,  the 
corresponding  Gray  Code  representation  for  3a  is  found  to  be  (0010)  ,  therefore  the 
recovered  data  state  will  be  Z?p- 

The  decoder  implementation  can  be  depicted  as  a  matrix  structure  in  a  manner 
similar  to  that  shown  for  the  transmitter  encoder  in  Chapter  II.  For  8-DPSK  and 
16-DPSK  modulated  signals,  the  decoder  matrices  are  given  in  Figure  3.2  and  Figure 
3.3  respectively.  For  example,  assuming  the  same  conditions  as  in  the  above 
paragraph,  by  using  S6  and  S3  as  a  column  and  row  entries  for  the  matrix  given  in 
Figure  3.3,  the  recovered  data  state  from  the  decoder  is  determined  as  D^  (0010). 
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Fiszure  3.2     Decoder  matrix  for  S-DPSK  sicnals. 


B.       COMPUTER  GENERATION  OF  RECEIVER  OUTPUTS 

Appendix  C  includes  the  results  of  the  computer  evaluation  of  the  M-ary  DPSK 
receiver.  In  addition  to  generating  the  decoder  matrix  in  decimal  form,  the  output  for 
this  program  includes  all  possible  combinations  of  the  decoder  outputs. 
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Figure  3.3    Decoder  matrix  for  16-DPSK  signals. 
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C.   RECEIVER  LOCAL  OSCILLATOR  PHASE  ERRORS 

The  introductory  chapter  summarizes  the  studies  done  on  M-PSK  DBDR 
structures,  assuming  that  the  receiver  local  oscillator  has  no  phase  error.  In  this 
research,  a  fixed  receiver  phase  error  has  been  assumed  in  the  performance  analyses 
conducted,  and  its  effect  on  receiver  BER  is  presented  in  Chapter  IV  and  Chapter  V. 
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IV.  PERFORMANCE  ANALYSES  OF  8-PSK  AND  8-DPSK  DBDR'S  WITH 
LOCAL  OSCILLATOR  PHASE  ERROR 


In  previous  chapters,  DBDR's  are  introduced  and  their  operational  properties  are 
explained.  In  this  chapter,  performance  analyses  for  8-PSK  and  8-DPSK  will  be 
carried  out  in  detail. 

A.       8-PSK    DBDR    PERFORMANCE    ANALYSIS    IN    THE    PRESENCE    OF 
RECEIVER  LOCAL   OSCILLATOR  PHASE  ERROR 

The  possible  transmitted  signals  which  are  represented  as  vectors  in  Figure  1.2. a  , 

can  be  mathematically  expressed  as 


Sj(t)  =  V2E~T;  sin  [  27ifQt  +  9.(t)  ]  ,        i=  1,2,3...,8  ,  0<  t  <Ts  (4.1) 

where 


Es=  j"     Si2(t)dt   ,  i=l,2,..,8  (4.2) 

0 

T    =  svmbol  duration 
s 

fQ  =  carrier  frequency  (Hz.) 

0.(t)  indicates  the  phase  of  the  transmitted  PSK  signal,  which  for  the  8-PSK  case 
becomes 


'  a  +  (i-1)  7C/4  i=  1,3,5,7 


Oi(t)  = 


0<t<T  (4.3) 


irt/4  -  a  i=  2,4,6,8 


The  angle  a   is   variable  as  shown  in  Figure  1.2. a;  however,  the  constraint 
0  <  a  <  7T/4  must  be  observed.   Equation  4.1  can  be  written  by  using  a  suitable 
orthonormal  set  of  functions  as 

s.(t)=  cp^t)  VES  cos9.(t)  +  <p2(t)  VES  sine.(t)  (4.4) 
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i=l,2 8   ;  0<t<T 


where  (p^t)  and  (p2(t)  are  the  orthonormal  function  set  which  by  applying   the 
Gram-Schmitt  orthonormalization  procedure 
[Ref.  3:  p.  198] 
can  be  shown  to  be 


q>j(t)  =  V2/T  sin  27tf0t 

(p2(t)  =  V2/Tcos  27tf0t  ,  0<  t  <Ts  (4.5) 

The  received  signal  depends  upon  the  channel  model  and  the  transmitted  signal 
as  specified  by  Equation  4.1.  In  this  research,  an  Additive  White  Gaussian  Noise 
(AWGN)  interference  model  has  been  assumed  so  that  the  received  signal  can  be 
written  as 

r(t)=s.(t)  +  n(t)  (4.6) 

where  n(t)  is  a  sample  function  of  a  white  Gaussian  noise  process  having  a  constant 
(one  sided)  power  spectral  density  function  equal  to  NQ  watts/Hz.  The  receiver 
processes  r(t)  over  an  interval  of  duration  Ts  seconds  and  a  decision  is  made  at  the  end 
of  the  interval.   This  procedure  continues  in  a  repetitive  fashion  as  time  progresses. 

Assuming  all  signals  are  equally  likely  to  be  transmitted,  the  probability  of  a  bit 
being  correct  for  M-PSK  signal  transmission  is  written  as 

M 

Pr{  bit  correct  }  =  l/M  £   Pr(  bit  correct  /  s.(t)  transmitted  }  (4.7) 

i=l 

The  calculations  for  the  conditional  probabilities  (Pr{  bit  correct  /  s.(t)  transmitted  }, 
i=  1,2,....M  )  have  been  determined  by  Myers  and  Bukofzer  [Ref.  7:  pp.  6-12],  for  the 
8-PSK  case  with  no  receiver  local  oscillator  phase  error.  In  the  sequel,  the  same 
analysis  methods  are  utilized  to  obtain  the  receiver  BER  in  the  presence  of  a  local 
oscillator  phase  error.  A  local  oscillator  phase  error  is  modeled  by  the  angle  z  in  the 
unmodulated  receiver  carriers  (Pj(t)  and  <p2(t)  of  Figure  1.6.  These  are  denoted  <p,'(t) 
and  (p2'(t),  namely 
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(pjXt)  =  V2/T  sin(  2icf0t  +  s) 

q>:'(t)  =  V2/T  cos  (  2:rf0t  +  €  ),         0  <  t  <  Ts  (4.8) 

Assuming  signal  Sj(t)  was  transmitted,  this  produces  outputs  rki  from  the  receiver's 
correlators,  where 

rk,i  =  1 1  tf«)  +  n^  1  %  w dt   •    i=  1'2-8  ■    k=  1-2  <4-9) 

o 

For  example  when  Sj(t)  is  transmitted,  from  Equation  4.4,  9.(t)  becomes  just  a  and 

r1  1  =  VES  cos  (a-e)  +  n, 

r2  j  =  VES  sin  (a-e)  +  n2  (4.10) 

where 


Ts 


nk  =  j"  n(t)  (pk'(t)  dt  ,        k=  1,2  (4.11) 

0 

Observing  that 

E  [n,  n2]  =  E  [  j"  n(t)  ^'(t)  dt  x     J  n(t)  <p2'(t)  dt  ]=  0  (4.12) 

0  0 

This  means  that  nl  and  n2  are  uncorrected  random  variables.  Since  n^  and  n2  are 
Gaussian  random  variables  (having  identical  variances)  [Ref.  7:  p. 7],  they  are  also 
statistically  independent. 

The  goal  of  the  receiver  of  Figure  1.6  is  to  correctly  recover  the  LB,  the  MB,  and 
the  RB,  transmitted  via  the  signal  s.(t).   By  introducing  the  more  compact  notation 

Pr{  bit  correct  /  s.(t) }  =  Pr{  bit  correct  /  s.(t)  transmitted  }   ,  i=  1,2.. ,8         (4.13) 
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Pr{  LB  correct  /  s.(t)  }  =     < 


the  probabilities  of  such  correct  bit  recovery  are  given  in  general  by 

i 

Pr{ru>0/^(t)}  i=  1,2,8,7 
Pr{rn<0/&(t)}  i=  3,4,5,6 
Pr{  r1}i2  -  r242>0  /  s.(t) }  i-  1,4,5,8 

I     Pr{rl,i2"r2,i2<°/Si(t))  i=   2<3'6'7 

Pr{r2i>0/s.(t)}  i=  1,2,3,4 
Pr{r2.<0/s.(t)}      i=  5,6,7,8 


Pr{  MB  correct  /  s.(t)  }  =   ( 


Pr{  RB  correct  /  Sj(t)  }  =    < 


From  Equations  4.8-4.16,  we  can  write 


(4.14) 


(4.15) 


(4  16) 


Pr{  LB  correct  /  s^t) }  =  Pr{  r1  x  >  0  /  Sj(t)  } 

=  Pr{  rij  >  -  VE^  cos(a-e)  /  s^t) } 
=  1  -  Q  [  V2  Es/N0  cos(a-e)  ] 

where 


oo 


Q(x)  =  J"    [  1/V2ir]  exp(-u2/2)  du 
x 


(4.17) 


(4.18) 


Similarly 


Pr{  RB  correct  /  Sj(t) }  =  Pr{  r2  j  >  0  /  Sj(t)  } 

=  Pr{  n2>  -  VE^  sin(a-e)  /  s1(t) } 
=  1  -  Q  [  V2  Es/N0   sin(a-e)  ] 


(4.19) 


In  order  to  evaluate  the  conditional  probability  of  the  MB  being  correct,  the  following 
steps  are  required.    From  Equation  4.15 


Pr{  MB  correct/  s.(t)  }    = 


Pr{rlfi2.r2.2>0/s.(t)} 
P'{(VV-(rl,i+V>0/si(t)} 


(4.20) 
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By  defining  variables  y,  .  and  v., .  as 

i.  .1  *   —  .1 

>U  =     rU  +  r2,i 

>*2.i  =     rl,i"r2,i  (4-21> 

then  from  Equation  4.9,  when  Sj(t)  is  transmitted, 

yl  {  =  V  Es  [  cos  (a-c)  +  sin  (a-e)  ]+  n{  +  n2  (4.22) 

y21  =  V  Es  [  cos  (a-e)  -  sin  (a-e)  ]+  nL  -  n2  (4.23) 

Clearly  y,  2  and  y2  j  are  Gaussian  random  variables  with  statistical  means  given  by 

E{  yl  {}  =  v^  [  cos  (a-e)  +  sin  (a-e)  ]  (4.24) 

E{  y2  i)  =  V~E^  [  cos  (a-e)  -  sin  (a-e)  ]  (4.25) 

and  variances  given  by 

Var{yM}  =  Var{  y2>1}  =  NQ  (4.26) 

By  observing  that 
E  (  &i,r  E I  ylflD  (y2,r  e  [  y2>1]) }  =   E  {  (n,  +  n2)  (nr  n2) } 

=    E  {  (ni)2-(n2)2  }  =  0  (4.27) 

where  nj  and  n2  are  defined  in  Equation  4.11.  Since  yl  j  and  y2  j  are  uncorrected 
Gaussian  random  variables,  they  are  also  statistically  independent.  Therefore, 
Equation  4.20  becomes 

Pr{  MB  correct  /  s^t) }  =    Pr{  yl  1   y2 1  >0  /  Sj(t) } 

-  Pr{  [  yM  >  0  /  s2(t)  ]  and  [  y2?1  >  0  /  s^t)  ]  }  + 

+  M  [  yM  <0  /  Sj(t)  ]  and  [  y2'i<0  /  Sl(t)  ]  }  + 

=  Pr{  y1(1>  0  /  Sl(t)  }  Pr{  y2>1  >  0  /  s^t)  }  + 

+  Pr{  yM  <  0  /  Sj(t)  }  Pr{  y2'1  <  0  /  s^t)  }  (4.28) 
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From  Equations  4.22  -  4.26,  we  have 


Pr{  yM >0  /  Sl(t)  }  =  Q  [  E{y14}  /  Vnq  ]  (4-29) 

Pr{  y2?1  >  0  /  Sl(t)  }  -  Q  [  E{yv}  /  V>^  1  <4-30) 


so  that  Equation  4.28  becomes 


Pr{  MB  correct  /  Sj(t) }  =  1  -  Q{  V  Es/N0  [  cos(a-e)  +  sin(a-e] } 

-  Q{  V  Es/N0  [  cos(a-c)  -  sin(a-e] } 

-  2  Q{  V  Es/N0  [  cos(a-e)  -  sin(a-e]  }  . 


.  Q{  V  Es/N0  [  cos(a-e)  +  sin(a-e] }  (4.31) 

Assuming  now  transmission  of  s2(t),  a  similar  analysis  procedure  is  performed 
but  now  Equations  4.14  -  4.16  become 

02(t)  =  7T/2  -  a  (4.32) 

rl  2  =   "^s  S*n  (a  +  £)  +  nl 

r22  =  VE^cos(a  +  £)  +  n2  .  (4.33) 

where  n,  and  n2  are  defined  by  Equation  4.11.   From  Equations  4.14  and  4.33 

Pr{  LB  correct  /  s2(t) }  =  Pr{  r{  2>0  /  s2(t) } 

=  Pr{  n2>  -  jEs  sin(a  +  £)  /  s2(t) } 

=  1  -  Q  [  V2  Es/N0   sin(a  +  e)  ]  (4.34) 

and 

Pr{  RB  correct  /  s2(t) }  =  Pr(  r2  2>0  /  s2(t) } 

=  Pr{  n2>  -  JTs  cos(a  +  e)  /  s2(t) } 

=  1  -  Q  [  V2  Es/N0  cos(a  +  e)  ]  (4.35) 

For  calculations  involving  the  MB,  one  additional  random  variable  needs  to  be  defined, 
namely 

y3li  ■  hi  -  hi      i=  l-2 8  <4-36) 
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Assuming  transmission  of  s,(t),  Equation  4.36  becomes 

-'3,2   =     r2,2  '  rl,2 

=    v'  Es  cos  (a  +  s)  +  n2  -  V  Es  sin  (a  +  s)-nj  (4.37) 

where 

E  (y3  2)  =  v^  cos  (a  +  £)  -  V~E^  sin  (a  +  e)  =  E  {y,  {}  (4.38) 

Since 

Pr{  MB  correct  /  s2(t) }  =    Pr{  (r{  {)2  -  (r2>2)2<  0  /  s2(t) }  (4.39) 

"   Pr^rl,'l-r2,2)'(rl,l   +   r2,2)<0    !*&)) 

by  using  Equations  4.21  and  4.36 

Pr{  MB  correct  /  s2(t)  }  =    Pr{  y{  2  y3  2<  0  /  s2(t)  }  (4.40) 

Furthermore 

E  (  (y1>2-  E  {  y1?2})  (y3)2-  E  {  y3>2}) }  =    E  {  (n{  +  n2)  (nr  n2) } 

=    E  {  (ni)2-(n2)2  }  =  0  (4.41) 

so  that  yt  2  and  y3  2  are  statistically  independent  Gaussian  random  variables,  and 
therefore  Equation  4.40  becomes 

Pr{  MB  correct  /  s2(t) }  =  Pr{  y3?2  y1?2  <  0  /  s2(t) } 

=  Pr{  [  yj>2>0  /  Sl(t)  ]  and  [  y3j2>0/  s2(t)  ]  }  + 

+  Pr{  [  y^2  <  0  /  s2(t)  ]  and  [  y3^2  <  0  /  s2(t)  ]  }  + 

=  Pr{  y1>2'> 0  /  s2(t)  }  Pr{  y3>2>  0  /  s2(t) }  + 

4-  Pr{  yx  2<0  /  s2(t)  }  Pr{  y3'2<0  /  s2(t)  }  (4.42) 

By  using  a  method  similar  to  that  used  in  the  derivation  of  Equations  4.29  and  4.30, 
for  Vj  2  and  y3  2  ,  Equation  4.42  becomes 
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Pr{  MB  correct  /  s2(t) }  -   1  -  Q{  V  ES/NQ  [  cos(a  +  s)  +  sin(a4-g]  }  - 

-  Q{  V  Es/N0  [  cos(a  +  s)  -  sin(a  +  e]  } 

-  2  Q{  V  Es/N0  [  cos(a  +  £)  +  sin(a  +  e]  }  . 

.  Q{  J  Es/N0  [  cos(a  +  e)  +  sin(a  +  £]  }  (4.43) 

Analyses  can  be  carried  out  for  probability  of  correct  reception  assuming  every 
possible  signal  transmission.  This  is  however  not  necessary  because  of  inherent 
symmetries  in  the  problem  .  In  conditions  where  symmetries  are  not  present,  previous 
methodologies  can  be  used  to  obtain  results  as  follows.  Assuming  transmission  of 
s3(t),  Equations  4.3  and  4.9  become 

93(t)  =  re/2  +  a  (4.44) 

rj  3  =  VE"S  cos  (a-£+re/2)  +  nx  (4.45) 

=  -  VES  sin  (a-e)  +  n1 
r23  =  VES  sin(a-£  +  7t/2)  +  n2  (4.46) 

=  -  VES  cos  (a-£)  +  n2 

so  that  from  Equation  4.14 

Pr{  LB  correct  /  s3(t)  }  =  Pr{  r{  3  <  0  /  s3(t)  } 

=  Pr{  Hj  >  -  VES  sin  (a-e)  /  s3(t) } 

=  1  -  Q[  V2  Es/NQ  sin(a-c)  ]  (4.47) 

and  from  Equation  4.16 

Pr(  RB  correct  /  s3(t)  }  =  Pr{  r2  3  <  0  /  s3(t) } 

=  Pr{  n2>  -  Ve"s  cos  (a-e)  /  s3(t) } 

-  1  -  Q[  V^ETNq  cos(a  +  e)  ]  (4.48) 

Assuming  transmission  of  s4(t),  Equations  4.3  and  4.4  become 


(4.49) 
(4.50) 


e4(t)  ■ 

=  7t-  a 

ri,4   = 

«/Es  COS  (71 

:-a-e)  + 

ni 

= 

-  JE   cos 

v     s 

(a  +  c)  + 

n 
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r24  =  vEs  sin(TT-a-e)  +  n,  (4.51) 

=  -  JWS  sin  (a  +  c)  +  n2 

so  that  from  Equation  4.14 

Pr{  LB  correct  ,  s4(t)  }  =  Pr{  r1>4  <  0  /  s4(t)  } 

=  Pr{  nx>  -  VES  sin  (o  +  c)  /  s4(t)  } 

-  1  -  Q[  V2  Es/N0  cos  (a  +  c)  ]  (4.52) 

and  from  Equation  4.16 

Pr{  RB  correct  /  s4(t)  }  =  Pr{  r2  4  >  0  /  s4(t)  } 

=  Pr{  n2>  -  VE"S  sin  (a  +  £)  /  s4(t) } 

-  1  -  Q[  «/2E/NQ  sin  (a  +  e)  ]  (4.53) 

Also  from  Equation  4.15 

Pr{  MB  correct  /  s4(t)  }  =    Pr{  (rI>4)2  -  (r2>4)2>  0  /  s4(t)  }  (4.54) 

=  Pr{  (rlf4  -  r2f4)  (r1>4  +  r24)>0  /  s4(t) } 

Using  Equations  4.21  and  4.36  yields 

>'l,4  =     '1,4  +   r2,4  (4-55) 

>'2,4  =     rl,4"r2,4  (4-56) 

and  using  a  method  similar  to  that  used  in  the  derivation  of  Equation  4.27  results  in 

Pr{  MB  correct  /  s4(t) }  =  Pr{  y1>4  y2>4  >0  /  s4(t) } 

=  Pr{  [  y'1?4>0  /  s4(t)  ]  and  [  y2j4>0  /  s4(t)  ]  }  + 

+  Pr{  [  y^4  <  0  /  s4(t)  ]  and  [  y2'4  <  0  /  s4(t)  ]  }  + 

=  Pr{  yM>  0  /  s4(t)  }  Pr{  y2>4>  0  /  s4(t) }  + 

+  Pr{  yM<0  /  s4(t)  }  Pr{  y2]4<0  /  s4(t)  }  (4.57) 

which  becomes 
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Pr{  MB  correct  /  s4(t) }  =  1  -  Q  [  *J2  ES,'NQ  [  cos(a  +  c)  -  sin(a  +  £)  ]  } 

-  Q  {  J  2  Es/N0  [  cos(a  +  £)  +  sin(a +  £)•.]}  + 
+  2  Q  {  V2  E./Nn  [  cos(a  +  e)  -  sin(a  +  s)  ] }  . 


.  Q  {  V2  Es/N0  [  cos(a  +  e)  +  sin(a  +  e)  ]}  (4.58) 

For  transmission  of  s5(t),  the  pertinent  equations  are 

05(t)  =  7t+  a  (4.59) 

r,  5  =  VES  cos  (7i  +  a-£)  +  rij  (4.60) 

=  -  «fEs  cos  (a-£)  +  n, 
r2  5  =  V^  sin(7t  +  a-£)  +  n2  (4.61) 

=  -  v^l  tin  (tt-fi)  +  n, 


so  that 

Pr{  MB  correct  /  s5(t)  }  =    Pr{  (r1>5)2  -  (r2>5)2>  0  /  s5(t)  }  (4.62) 

asPrHr1>5-r2>5)'(r1}5  +  r2>5)>0  /  s5(t) } 

which  becomes 

Pr{  MB  correct  /  s5(t)  }  =  Pr{  y1>5  y2>5  >0  /  s5(t)  } 

-  Pr{  [  y'1>5>0  /  s5(t)  ]  and  [  y2?5>0  /  s5(t)  ]  }  + 
+  Pr{  [  yM<0  /  s5(t)  ]  and  [  y2^5<0  /  s5(t)  ]  }  + 
=  Pr{  y1?5'>  0  /  s5(t) }  Pr{  y^>  0  /  s5(t)  }  + 

+  Pr(  y^5<0  /  s5(t) }  Pr{  y2'5<0  /  s5(t) }  (4.57) 

by  observing   that  y.  5  and  y2  5  are  statistically  independent  random  variables  making 
it  possible  to  apply  the  similar  procedure  shown  above.   The  result  becomes 

Pr{  MB  correct  /  s5(t)  }  =  1  -  Q  {  J 2  ES/NQ  [  cos(a-E)  +  sin(a-c)  ]  }  - 

-  Q  {  V2  Es/N0  [  cos(a-£)  -  sin(a-e)  ]}  + 

+  2  Q  (  V2  Es/N0  [  cos(a-£)  +  sin(a-£)  ]  }  . 

.  Q  {  V2  Es/N0  [  cos(a-£)  -  sin(a-£)  ]}  (4.64) 
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As   indicated  before,   because  of  the  inherent  symmetry  in  the  problem,   the 
following  equations  hold  for  conditional  LB  correct  probabilities  (see  Figure  1.2.b  ) 

Pr{  LB  correct  /  s^t)  }  =  Pr{  LB  correct  /  s5(t)  }  (4.65) 

Pr{  LB  correct  /  s2(t)  }  =  Pr{  LB  correct  /  s6(t)  }  (4.66) 

Pr{  LB  correct  /  s3(t)  }  =  Pr{  LB  correct  /  s?(t)  }  (4.67) 

Pr{  LB  correct     s4(t) }  =  Pr{  LB  correct  /  sg(t)  }  (4.68) 

and  similarly  for  the  RB 

Pr{  RB  correct     Sj(t)  }  =  Pr{  RB  correct  /  s5(t)  }  (4.69) 

Pr{  RB  correct  /  s2(t)  }  =  Pr{  RB  correct  /  s6(t)  }  (4.70) 

Pr{  RB  correct  /  s3(t)  }  =  Pr{  RB  correct  /  s?(t)  }  (4.71) 

Pr{  RB  correct  /  s4(t) }  =  Pr{  RB  correct  /  sg(t)  }  (4.72) 

and  finally  for  the  MB 

Pr{  MB  correct  /  Sj(t)  }  =  Pr{  MB  correct  /  s3(t)  }  (4.75) 

Pr{  MB  correct  /  s2(t)  }  =  Pr{  MB  correct  /  sg(t)  }  (4.76) 

Pr{  MB  correct  /  s4(t)  }  =  Pr{  MB  correct  /  s6(t)  }  (4.77) 

Pr(  MB  correct  /  s5(t)  }  =  Pr{  MB  correct  /  s?(t)  }  (4.78) 

By  using  now  Equation  4.7  for  the  unconditional  probability  of  the  LB,  MB  and 
RB  being  correct,  the  above  results  yield 

Pr{  RB  correct  }  =  Pr{  LB  correct  } 


=  1  -  1/4  [  Q[  V2E7N0  cos(a-e)  ]  + 

+  Q[V2Es/N0sin(a  +  £)]  + 

+  Q[  V2ES/N0cos(a  +  E)] 

+  Q[  J2EJWQ  sin(a-e)  ]  ]  (4.79) 


and 


Pr{  MB  correct }  =  1  -  1/2  [  Q  {  J7  Es/NQ  [  cos(a-£)  +  sin(a-e)  ]  }  + 
+  Q  {  V2  Es/N0  [  cos(a  +  e)  +  sin(a  +  £)]}  + 
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+  Q  {  V2  Es/N0  [  cos(a-e)  -  sin(a-e)  ] }  + 
+■  Q  {  V2  Es/N0  [  cos(a  +  £)  -  sin(a  +  e)  ]}  ]+ 
+  Q  {  J 2  Es/N0  [  cos(a-e)  +  sin(a-e)  ] }  . 
.  Q  {  V2  Es;N0  [  cos(a  +  e)  +  sin(a  +  c)  ]}  + 
+  Q  {  V2  Es/NQ  [  cos(a-c)  -  sin(a-e)  ]  }  . 
.  Q  {  V2  Es/N0  [  cos(a  +  e)  -  sin(a  +  €)]}]  (4.80) 

For  8-PSK,the  Signal  to  Noise  Ratio  (SNR)  has  been  defined  by  Equation  1.6,  so 
that  using  Equations  4.79  -  4.80,  the  BER  denoted  PE  can  be  expressed  as  a  function 
of  SNR,  a  and  £  as  follows 

PE=  1/3  [  Pr{  LB  correct  }  +  Pr{  MB  correct  }  +Pr{  RB  correct  }  ] 


=  1  -  1/6  {  Q[  2  VSNR  cos(a-£)  ]  +  Q[  2  VSNR  sin(a  +  e)  ]  + 


+    Q[  2  VSNR  cos(a  +  £)  ]  +  Q[  2  VSNR  sin(a-£)  ]  + 


+  Q  {  2  V  SNR  [  cos(a-£)  +  sin(a-c)  ]  }  + 
+    Q  {  2  V  SNR  [  cos(a  +  £)  +  sin(a  +  £)]}  + 
+  Q  {  2  V  SNR  [  cos(a-£)  -  sin(a-£)  ]  }  + 
+  Q  {  2  V  SNR  [  cos(a  +  £)  -  sin(a  +  £)]  }  }  + 
+  Q  {  2  J  SNR  [  cos(a-£)  +  sin(a-£)  ]  }. 
.   Q  {  2  V  SNR  [  cos(a  +  £)  +  sin(a  +  £)]  }  + 
+  Q  {  2  V^NR  [  cos(a-£)  -  sin(a-E)  ] }. 


.  Q  {  2  V  SNR  [  cos(a  +  £)  -  sin(a  +  £)  ]  }  (4.81) 

From  Equation  4.81,  the  PE  can  be  computed  for  various  values  of  SNR,  a,  and 
£,  with  the  results  used  to  optimize  the  DBDR  performance.  For  £=0°  (i.  e.  ,  no 
receiver  local  oscillator  phase  error),  values  of  PE  vs  SNR  have  been  plotted  (see 
Figure  4.1)  showing  complete  agreement  with  the  results  obtained  by  Myers  and 
Bukofzer  [Ref  7:  p.  11]  Figure  4.1  includes  also  the  PE  for  various  values  of  £,  as  a 
function  of  SNR. 

B.       8-DPSK    DBDR    PERFORMANCE    ANALYSIS    IN    THE    PRESENCE    OF 
RECEIVER  LOCAL   OSCILLATOR  PHASE  ERROR 

The  advantages  and  disadvantages  of  differential  encoding  and  decoding  have 

been  previously  explained  in  the  introductory  chapter.    In  this  section,  performance  of 

the  8-DPSK  DBDR  is  presented.   As  explained  earlier  the  8-DPSK  decoder  utilizes  the 
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Figure  4. 1     PE  Versus  SNR  of  8-PSK  DBDR  for  Various  Values  of  £  (a  =  22.5°). 
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LB.  MB  and  RB.  which  are  recovered  via  direct  bit  detection  methods,  and  performs  an 
inverse  mapping  from  which  the  original  data  stream  is  delivered  to  the  intended  user 

Before  proceeding  with  the  performance  calculations,  the  following  notation  will 
be  introduced.   The  probabilities  of  a  bit  being  correct  or  being  in  error  will  be  denoted 

Ck  =  Pr{  k  correct }    ,        Ek  =  Pr{  k  error  }  (4.82) 

respectively,  where  k=  1  for  LB,  k=  2  for  MB,  k=  3  for  RB.  By  introducing  the 
notation  given  by  Equation  4.82,  the  analysis  for  the  performance  of  the  8-DPSK 
DBDR  can  be  expressed  in  more  compact  form.  For  example  E2  means  the 
probability  that  the  MB  is  in  error. 

In   order  to    obtain  the  probability   of  a  bit   being   correct  (with  differential 
decoding)  we  must  fust  evaluate 

Pr{  RB  correct }  =  1/2  Pr{  RB  correct  /  RB=  1  }  + 

+    1/2  Pr{  RB  correct  /  RB=  0  }  (4.83) 

where  by  further  conditioning,  we  equivalently  have 

1        1 
Pr{  RB  correct/ RB  =  m  }  =  £      £        Pr{  RB  correct  /  RB  -  m,MB  =  i,LB  =  j  } . 

i  =  0    j  =  0 

.  Pr{  MB  =  i,  LB  =  j  }         m  =  0,l  (4.84) 

Due  to  the  assumption  of  equally  likely  probabilities  of  having  O's  and  l's  in  the  data, 

Pr{  MB  =  i,  LB  =  j  }  =  1/4  i,j=  0,1  (4.85) 

so  that  Equation  4.83  becomes 

1      1 
Pr{  RB  correct }  =  1/8  £    £  Pr{  RB  correct  /  RB  =  0,MB  =  i,LB  =  j  }  + 

i  =  0  j  =  0 
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1     1 

+  1/8  y     £  Pr{  RB  correct  /  RB=  l,MB  =  i,LB  =  j  }  (4.86) 

i  =  0    j  =  0 

The  DBDRs  recover  from  the  received  signal,  each  bit  in  the  data  state  (for  PSK 
modulation)  or  each  bit  in  the  symbol  state  (for  DPSK  modulation)  independently. 
This  property  of  the  DBDR  is  an  advantage  for  BER  analyses.  Because  certain 
consecutive  bit  errors  in  the  past  and  presently  recovered  symbol  states  do  not 
necessarily  result  in  errors  in  the  recovered  data  states,  these  bit  errors  in  successive 
symbol  states  must  be  considered.  That  is,  even  though  symbol  errors  may  occur,  the 
delivered  data  may  stillbe  correct.  Therefore,  these  type  errors  have  to  be  accounted 
for  as  resulting  in  a  correct  reception  and  considered  in  the  computation  of  the 
probability  of  a  bit  being  correct.  Figure  4.2  shows  one  example  of  the  situation  being 
discussed.  Here  the  data  state  (0,0,0)  is  assumed  to  occur  and  the  symbol  state  (0,0,0) 
is  assumed  to  have  been  previously  transmitted.  This  means  that  according  to  the 
encoder  matrix  shown  in  Chapter  II,  the  symbol  state  (0,0,0)  must  now  be  generated 
for  transmission.  Figure  4.2  further  shows  the  receptions  of  the  symbol  states  which 
even  though,  most  contain  errors,  produce  recovered  data  states  in  which  the  RB  is 
correctly  delivered.  Observe  that  there  are  32  such  combinations  which  yield  a  correct 
0  for  the  RB.  Therefore  notation  has  to  be  introduced  in  order  to  account  for  error 
patterns  of  the  above  type,  namely  Pr{  RB  correct,  (k,l,m)  /  RB  =  0,MB  =  i,LB  =  j  }, 
where  k,l,m=  0,1.  The  meaning  of  this  expression  is  the  probability  of  the  RB 
correctly  recovered  given  that  the  data  state  was  (j,i,0)  and  the  delivered  data  state  is 
(k,l,m)  (where  m=0  in  this  case,  otherwise  the  RB  could  not  be  correct  due  to  the 
transmission  of  the  RB  =  0).   Therefore 

Pr{  RB  correct  /  RB  =  0,MB  =  0,LB  =  0  }  = 
1       1 
=  X      Z   Pr(  RB  correct,(k,l,0)  /  RB  =  0,MB  =  0,LB  =  0  }  (4.87) 

k  =  0  1  =  0 

Since  each  element  of  the  symbol  state  is  independently  recovered,  we  see  from  Figure 
4.2  that,  for  example 
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CORRECT  RECEPTION  COMBINATIONS  IN  PRESENCE  OF  ERROR 
^MHHMHhHHHHfiHHHHHHHHHHHHHHHHHMHHHHHHhlHHHHHMHHMHHHHHmiHHHHHHHHHHHHhlMm 


INPUT 
0     0 


PREVIOUS  PRESENT  PREVIOUS  PRESENT  OUTPUT 

TRANSMISSION  TRANSMISSION       RECEPTION       RECEPTION 


Figure  4.2     Example  of  Error  Patterns  for  8-DPSK  System  Where  Data  State. 
Previous  Transmission,  and  Present  Transmission  are  000  with  RB  Correctly  Recovered. 
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Pr{  RB  correct,(0,0,0)  /  RB  =  0,MB  =  OXB  =  0  }  =  CI2  C22  C32  + 
+  CI2  C22  E32  +  CI2  E22  C32  +  CI2  E22  E32  +  El2  C22  C32  + 

+  El2  C22  E32+  El2  E22  C32+  El2  E22  E32  (4.88) 

From  Equations    4.79   and  4.82  -   4.86,   it  is   clear  that  Cl=    C3   and  El=    E3. 
Therefore,  simplification  of  Equation  4.88  is  possible,  thus  yielding 

Pr{  RB  correct,(0,0,0)  /  RB  =  0,MB  =  0,LB  =  0  }  =  CI4  C22+  2  CI2  C22  El2  + 

+  CI4  E22+  2  CI2  El2  E22+  El4  C22+  El4  E22  (4.89) 

A  similar  expression  must  be  obtained  for  other  combinations  of  values  of  k  and  1  as 
well  as  the  condition  that  the  RB  is  1.  This  produces  a  total  of  64  possibilities  of 
received  symbol  states  that  result  in  a  correct  RB.  A  computer  generated  table  was 
used  ,in  order  to  identify  the  equivalent  error  patterns  for  other  cases  involving  MB=  1 
and  LB=  1.  The  appropriate  listing  are  shown  in  Appendix  C.  Similar  analyses  must 
be  carried  out  in  order  to  determine  the  unconditional  probabilities  Pr{MB  correct}  as 
well  as  Pr{LB  correct}.    Such  analyses  demonstrate  that  due  to  the  inherent  symmetries 

Pr{  LB  correct  }  =  Pr{  MB  correct  }  =  Pr{  RB  correct  }  ■  PE  (4.90) 

so  that  using  Equation  4.86  along  with  similar  forms  obtained  for  k=  1  and  1=  1  yield 

PE  =  1  -  1/2  [  2C14  C22  +  9  CI2  C22  El2+  CI4  E22  +  3  CI4  E2  C2  + 
+  10  CI2  El2  E22+  2  El4  C22  4-  El4  E22  +  17  CI2  C2  El2  E2  + 
+  CI3  C22  El  +  3  CI3  El  E22  +  2  CI  El3  C22  +  2  CI  El3  E22  + 
+  3  C2  El4  E2    +  4  CI  C2  El3  E2  +  4  CI3  C2  El  E2  ]  (4.91) 

In  Figures  4.4  and  4.5,  PE  (given  by  Equation  4.91  )  vs  SNR  have  been  plotted  for 
various  values  of  a  and  £  respectively. 

C.  NUMERICAL  RESULTS  AND  DISCUSSION  OF  RESULTS  FOR  8-PSK 
AND  8-DPSK  DBDR'S  WITH  RECEIVER  LOCAL  OSCILLATOR  PHASE 
ERROR. 

Previous  sections  of  this  chapter  have  focused  on  the  performance  of  8-PSK  and 

8-DPSK  DBDRs  operating  in  the  presence  of  a  receiver  local  oscillator  phase  error. 

Without  any  phase  error,  we  have  shown  the  performance  plots  of  8-PSK  and  8-DPSK 
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Figure  4.3     PE  vs  SNR  Comparison  Plot  for  8-PSK  and  S-DPSK  DBDR 

(a=  22.5°  ,  e=  0°). 
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DBDRs  on  Figure  4.3,  thus  demonstrating  the  performance  penalty  of  the  latter 
receiver.  While  usage  of  8-DPSK  provides  good  bandwidth  efficiency  (same  as  8-PSK) 
and  relatively  good  noise  immunity,  the  BER  of  S-PSK  is  lower  than  that  of  8-DPSK 
for  £=  0°,  as  expected 

The  value  of  using  8-DPSK  versus  8-PSK  can  be  appreciated  more  fully  when 
comparing  the  performance  plots  of  Figures  4.1  and  4.5.  These  show  that  if  the 
receiver  has  a  local  oscillator  phase  error  restricted  to  0  ^  e  ^  66.5°,  for  similar  SNR 
values,  the  BER  of  the  8-DPSK  receiver  is  lower  than  that  of  the  8-PSK  receiver. 
Since  such  phase  errors  can  be  typically  expected  to  be  small,  this  is  a  significant 
advantage.  However  it  must  be  pointed  out  that  for  phase  errors  of  more  than  10°, 
both  receivers  exhibit  severely  degraded  performance  that  would  make  both  systems 
unsuitable  for  operation  in  a  practical  system.  The  8-DPSK  receiver  does  not  exhibit 
complete  insensitivity  to  phase  errors  as  encountered  in  phase  measurement  receivers. 
Only  when  a  complete  phase  revesal  takes  place  (i.e.,  c=  180°)  is  the  performance  of 
the  8-DPSK  DBDR  unaffected,  while  the  8-PSK  receiver  produces  a  nearly  useless 
output. 

The  angle  a  set  to  22.5°  is  optimum  in  minimum  BER  sense  for  the  8-PSK 
DBDR.  The  variations  in  PE  versus  a  are  shown  in  Figure  4.6,  for  e=  0°  for  8-PSK, 
and  a  similar  plot  of  PE  versus  a  is  plotted  for  various  values  of  SNR  and  e=  0°  on 
Figure  4.7  for  the  8-DSPK  case. 
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V.  PERFORMANCE  ANALYSES  OF  16-PSK  AND  16-DPSK  DBDR'S 
WITH  LOCAL  OSCILLATOR  PHASE  ERROR 


In  this  chapter,  DBDRs  will  be  analyzed  for  16-PSK  and  16-DPSK 
communication  systems  operating  in  the  presence  of  a  receiver  local  oscillator  phase 
error.  The  analysis  procedure  will  be  similar  to  that  carried  out  in  the  previous  chapter 
for  8-PSK  and  8-DPSK  systems. 

The  performance  of  16-PSK  DBDR's  has  been  investigated  and  the  methodology 
for  carrying  out  the  performance  analyses  have  been  well-documented  by  Bukofzer 
[Ref.  8],  In  this  research,  these  analyses  will  be  extended  to  the  16-DPSK  system  as 
well  as  the  16-PSK  system  modified  to  operate  in  the  presence  of  a  receiver  local 
oscillator  phase  error. 

A.       16-PSK    DBDR    PERFORMANCE    ANALYSIS    IN    THE    PRESENCE    OF 
RECEIVER  LOCAL  OSCILLATOR  PHASE  ERROR 

The  possible  transmitted  signals,  which  are  represented  as  vectors  in  Figure  1.2.b, 

are  mathematically  given  by  Equation  4.1    with  1=    1,2, 3,... ,16.    The  same  notation 

introduced  in  Chapter  IV  will  be  used  to  represent  the  signals  using  the  orthonormal 

set  given  by  Equation  4.4.   Therefore  the  transmitted  signals  are  expressed  as 

sj(t)=  (pL(t)  VFcosa(t)  +  <p2(t)  VE;  sinGjCt)  (5.2) 

where  i=  1,2,..., 16,  0<t<Tg   and 

fa  +  (i-l)7t/8  i=  1,5,9,13 

P  +  (i-2)7i/8  i=  2,6,10,14 

9.(t)  =  {  (i+l)7l/8-P  1=3,7,11,15  (5.3) 

i  7i/8 -a  1  =  4,8,12,16 

with  a  and  P  as  shown  on  Figure  1.2.b 

Assuming  all  signals  are  equally  likely  to  be  transmitted,  for  the  16-PSK  case 
Equation  4.7  can  be  written  as 
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16 
Pr{  bit  correct  }  =  1/16  £    Pr{  bit  correct  /  s.(t)  transmitted  } 

i=l 


(5.4) 


The  notation  to  denote  the  four  bits  that  make  up  a  data  state  has  been  introduced  in 
Chapter  I.   Also  the  notation  introduced  by  Equation  4.13  is  used  here,  namely 

Pr{  bit  correct  /  Sj(t) }  =  Pr{  bit  correct  /  Sj(t)  transmitted  }   i=  1,2.., 16  (5.5) 

From  Figure  1.2.b,  the  conditional  probabilities  for  the  correct  reception  of  the  LB, 
MBL,  MBR.  RB  are  given  by 


Pr{  LB  correct  /  Sj(t)  }    =  < 


Pr{ru>0/si(t)}         1=1,2,3,4,13,14,15,16 

Pr{ri;<0/s.(t)}         1=5,6,7,8,9,10,11,12 
y,  1,1  i 


(5.6) 


Pr{  r1?i2-r2ji2>0  /  s.(t)  }        i=  1,2,7,8,9,10,15,16 
Pr{  MBL  correct  /  Sj(t)  }  =i  (5.7) 

Pr{  r,  :2-r2i2<0  /  s.(t)  }         1=3,4,5,6,11,12,13,14 


Pr{  MBR  correct  /  s.(t)  }  =  < 


Pr{  RB  correct  /  s.(t)  }    =    < 


Pr{  VTh>  0  /  Sj(t)  }  i=  1,4,5,8,9,12,13,16 

Pr{  VTh<0  /  Sj(t)  }  1=2,3,6,7,10,11,14,15 

Pr{r24>0/s.(t)}         1=1,2,3,4,5,6,7,8 


(5.8) 


(5.9) 


Pr{  r7.<0/  s.(t)}         i=  9,10,11, 12,13,14,15,16 


where  the  r.  .  are  given  by  Equation  4.9  and 


2  ,     2  \2 


Th  =  (  V  V  )  -  4  (ri,ir2,i)    .        i-W ,16 


(5.10) 


The  analyses  methods  which  have  been  given  in  Chapter  IV,  can  be  applied  to 
determining  the  conditional  error  probabilities  for  LB,  MBL,  RB.  For  MBR,  the 
analysis  procedure  which  has  been  used  by  Bukofzer  [Ref.  8],  will  be  utilized.    In  order 
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to  define  the  conditional  probability  sets  of  Equation  5.8,  additional  steps  must  be 
completed  as  demonstrated  in  the  sequel.  If  s.(t)  transmitted,  the  MBR  is  correctly 
recovered  if  for  i=  1,4,5,8,9,12,13,16,  VTh>0  and  if  for  i=  2,3,6,7,10,1 1,13,14,  VTh<0. 
Even  though  the  r  .  are  conditional  Gaussian  random  variables,  because  of  the  highly 
nonlinear  form  of  Equation  5.10,  a  more  productive  approach  to  the  determination  of 
the  conditional  p.d.f.  of  Vy^  has  been  introduced  (see  Reference  8  for  details)  rather 
than  using  transformations  of  Gaussian  random  variables.  By  introducing  the 
following  notation,  namely 

rj  .=»  component  of  r(t)  along  <p,(t)  =  t  cos  x\  (5.11) 

r2j=  component  ofr(t)  along  (p2(t)  =  €  sin  t]  (5.12) 

it  can  be  demonstrated  that 

VTh  =  C4  cos  4ii  (5-13) 

where 
e  =  length  of  the  vector  =  V(r,  .)2  +  (r,  .)2  (5.14) 

r\  =  vector  angle  that  with  respect  to  the  positive  x  axis  (5.15) 

=  arctan  (r2f|  /  r^) 

Since  t4>  0,  it  is  clear  that  the  test  involving  Vj^  being  greater  than  or  less  than  zero 
is  now  equivalent  to  the  test  involving  cos4ti  in  the  same  manner.  This  implies  that 
the  p.d.f.  of  r\  must  be  obtained  which,  as  shown  in  Reference  8,  leads  to  integrals  that 
cannot  be  evaluated  in  closed  form  and  therefore  must  be  evaluated  via  numerical 
integration. 

By  utilizing  the  analysis  procedure  from  Reference  8  and  a  methodology 
introduced  in  the  previous  chapter  for  the  8-PSK  case,  the  conditional  probabilities  for 
the  correct  reception  of  a  bit  in  the  presence  of  receiver  local  oscillator  phase  error  is 
given  in  shorthand  form  by  the  following  expressions 


Pr{  LB  correct  /  Sj(t)  }  =  1  -  Q{  Jl  ES/NQ  f  (a.p.e)  },    i-  1,2,3... 16  (5.16) 

where 
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f  (a.p,£)  = 


cos  (a-e) 

i=l,9 

cos  (P-£) 

i=2,10 

sin(p  +  £) 

i=3,ll 

/   sin(a  +  £) 

i— 4,12 

sin  (p-£) 

i=5,13 

sin  (a-£) 

i=6,14 

cos  (p  +  £) 

i=7,15 

cos  (a  +  e) 

i=8,16 

(5.17) 


and  a  and  P  are  shown  on  Figure  1.2.b.   Furthermore 


Pr{  MBL  correct  /  s.(t) }  =  Q{-V  ES/NQ  (cos  y  +  sin  y)  .  Q(-VEs/NQ  (cos  y-sin  y) }  + 

+  Q{VES/N0  (cos  y-sin  y) }.  Q{VES/N0  (cos  y  +  sin  y) }     (5.18) 


where 


y  =   f(a-£)     i=  1,5,9,13 
\(p-£)     i=  2,6,10,14 


(5.19) 


or 


Pr{  MBL  correct  /  Sj(t) }  =  1  -  Q{-VES/NQ  (cos  y  +  sin  y) .  Q{-VES/NQ  (cos  y-sin  y) }  + 

4-  Q{ VES/NQ  (cos  y  +  sin  y)}  .   Q{ VES/NQ  (cos  y-sin  y)}    (5.20) 


where 


y  =   f(P-£)     i=  3,7,11,15 
l(a-£)     i=  4,8,12,16 


(5.21) 


For  the  calculations  involving  the  conditional  probability  for  the  MBR,  the  joint 
p.d.f.  of  the  random  variables  r\  and  I  has  been  determined  by  Bukofzer  [Ref.  8] 
namely 


fftl,.  (L  ,H)  =    |L|  fr     r  (L  cosH  ,  L  sinH)  + 

Ml  bi  .  ri,i'r2,i/s: 


(5.22) 


+  |L|  L     ,'1      (-L  cosH  ,  -L  sinH) 

rl,i'r2,i/Sj 
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which  must  be  integrated  to  produce  the  marginal  p.d.f.  of  T|.   This  results  in 

fy    (H  /  s.)  =  "l/2  exp(-Es;N0)  +  {  1  +    V  4Es/JrNQ  .  (5.23) 

'  \  {  exp[  ES/NQ   cos2  (H-y]   cos(H-^)  .  Q[  -Jl  Es/N0  (cosfH-fc.)  ]   } 

where 

^.  =  e.(t)-e  i=l,2....16  (5.24) 

where  0.(t)  is  given  by  Equation  5.3  Therefore,  the  probability  that  Vj^>0  involves 
computing  the  probability  that  cos4ri>0.  Observe  that  cos4r(>0  for  -7T/8  <  x\  <  Tt/8, 
3n/8 <x\< 57C/8,  7tt/8 <r\< 971/8,  and  1  Itt/8 <r\<  13tc/8.  Therefore  the  calculation  of 
the  conditional  probability  of  the  MBR  being  correct  given  that  s^t)  was  transmitted 
can  be  calculated  by  integrating  the  conditional  p.d.f  of  r\  over  the  appropriate 
regions,  namely 

Pr{  MBR  correct  /  Sj(t)  }  =  Pr{  cos4q>0  /  s1(t)  }  = 

7T/8  5TC/8  971/8 

-    J  fn/sx  tH/sJ  dH  +  J  fn/Si  [H/Sl]  dH  4-  J  fn/Si   [H/Sl]  dH  + 

-7C/8  3TC/8  771/8 

137C/8 
+  J  fn/Si  [H/sJ  dH 

llTt/8 

71/ 8 

=    J  {  fy    [H  /  sj  +  fn/Si  [(H  +  w/2)  /  sx]  +  fy     [(H  +  tc)  /  Sj]  +  (5.25) 

-7T/8 

+  fn/s   [(H  +  371/2)  /sJJdH 

By  utilizing  the  appropriate  form  for  the  conditional  p.d.f.  of  r\  when  s.(t)  is  assumed 
transmitted,  results  in 

Pr{  MBR  correct  /  s.(t)  }  =  1/2  exp(-Es/NQ)  + 

TT/8 

+  V  Es/rcN0   J    (  exp[-  Es/N0  sin2(H-X)]  cos(H-X)  .  {Q[-J2  ES/NQ  (cos(H-X)]  - 
-Tt/8 

64 


-  Q[v'2  Es/N0  (cos(H  +  X)]  }  +  (  exp[-Es/NQ  cos2(H-),)]  sin(H-X)  Q[V2  EJWQ    (sin(H-X)] 

-  Q[v2  Es,-N0  (sin(H  +  X)  ]  }  }  dH  (5.26) 

where 


X  = 


< 


fa  +  £  i  =  4,8,12,16 

a-s  i=  1,5,9,13  (5.27) 

P-e  i=  2,6,10,14 

P  +  e  i=3,7,ll,15 


Finally,  it  can  be  seen  that 

Pr{  RB  correct  /  s^t)  }  =  Pr{  LB  correct  /  s.(t)  }  (5.28) 

so  that  by  using  Equations  5.4  and  5.16-5.28  for  correct  reception  probabilities  of  the 
individual  bits,  the  overall  BER  can  be  obtained  by  using 

Pr{  bit  in  error  }  =  1  -  Pr{  bit  correct  }  (5.29) 

This  results  in  the  (unconditional)  probabilities 

Pr{  LB  in  error  }  =  Pr{  RB  in  error  } 


=  1/8  {  Q[  V2  E  'Nn  cos(a-c)]  +  Q[  J2  E5/Nn  cos(P-c)]  + 


+  Q[  V2  ES.'N0  sin(p  +  c)]  +  QU/2  Es/N0   sin(a  +  £)]  + 

+  Q[  V2  Es/NQ  sin(P-£)]  +  QU/2  ES/NQ   cos(P  +  £)]  + 

+  Q  [V2  E~/N0  cos(a  +  £)]  +    Q  [V2  E~/NQ   sin(a  +  £)]]  (5.30) 

and 

Pr{  MBL  in  error  }  =  1/2  +  1/4  [  Q  {-VE'Nn  [cos(a-c)  +  sin(a-c)]}  + 


+  Q  {-  VE~X0  [cos(a  +  £)  -  sin(a  +  £)]}  +  Q  {-JEJNq  [cos(P-£)  +  sin(P-£)]}  + 
+  Q  {-VE~\  [cos(P-£)  -  sin(P-£)]}  -  Q  {-JEJN0  [cos(P  +  £)  +  sin(P  +  £)]}  - 
-  Q  {-JE/SQ  [cos(P  +  £)  -  sin(p  +  £)]}  -  Q  {-VES/N0  [cos(a  +  £)  +  sin(a  +  £)]}  - 
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-  Q  (VE/Nq  [cos(a  +  £)  -  sin(a  +  e)]}  ]  +  1/2  [  Q  {-*/EjRQ  [cos(a  +  e)  +  sin(a  +  £)]} 
.  Q  {-VE~N0  [cos(a-e)  -  sin(a-e)]}  +  Q  {-Ve7N0  [cos(P  +  £)  +    sin(p  +  £)]}  . 
.  Q  {-v/^/N0  [cos(P  +  £)  -  sin(p  +  £)]}  -  Q  {-VE/NQ  [cos(a-e)  +  sin(a-e)]}  . 
.  Q  [VE^Nq  [cos(a-E)  -  sin(a-£)]}  -  Q  {-VE^o  [C0S(P"e)  +    sin(P-£)]}  • 


Q{-VElN:[cos(P-e).sin(P-€)]}] 


(5.31) 


For  the  MBR  the  result  appears  in  integral  form.   That  is 


71/ 8 


Pr{  MBR  in  error  }  =  1  -  1/2  V  Es/7iN0  J  {  exp{-Es/NQ  sin2[H-(a-c)]}  cos[H-(a-e)]. 

-71/8 

{  1/2  -  Q[-VE7N0  cos(H-(a-£))]}  +  exp{-Es/N0  cos2fH.  <a-t)]}  sin[H-(a-e)] 
1/2  -  Q[-VE7^  sin(H-(a-£))]}  +  exp{-Es/NQ  sin2[H-{P-£)]}  cos[H-(p-e)j 
1/2  -  Qt-VETN'o  cos(H-(P-£))]}  +  exp{-Es/NQ  cos2[H-(P-£)]}  sin[H-(p-£)] 
1/2  -  Q[-Ve7n^  sin(H-(p-£))]}  -  exp{-Es/N0  sin2[H  +  (P  +  £)]}  cos[H  +  (p  +  £)] 
1/2  -  Q[-yE7N"0  cos(H  +  (p +  £))]}  -  exp{-Es/N0  cos2[H  +  (P  +  £)]}  sin[H  +  (P  +  £)] 
1/2  -  Q[-VE7N^  sin(H  + (a +  £))]}  +  exp{-Es/NQ  sin2[H  +  (a  +  £)]}  cos[H  +  (a  +  £)] 
1/2  -  Q[-VE 'Nn  cos(H  +  (a  +  £)]}  +  exp(-E  /Nn  cos2[H  +  (a  +  £)]}  sin[H  +  (a  +  £)] 


1/2  -  Q[-Ve7n:  sin(H  + (a +  £))]}  }  dH 


(5.32) 


The  overall  BER  can  now  be  determined  from  Equations  5.30,5.31  and  5.32  by 
evaluating  the  final  result  for  PE,  the  bit  error  probability,  namely 


PE  =  1/4  [  Pr{  LB  in  error  }  +  Pr{  MBL  in  error  } 
+  Pr{  MBR  in  error  }  +  Pr{  RB  in  error  }  } 


(5.33) 


It  can  be  seen  that  Equation  5.33  consists  of  many  terms  that  by  themselves  do  not 
yield  as  much  insight  as  the  computed  values  of  PE  presented  in  plotted  form. 
Therefore  plots  of  BER  vs.  SNR  for  various  values  of  a  and  £,  are  given  as  Figures  5.1 
and  5.2.  While  a  and  P  can  be  set  independently,  from  an  optimization  point  of  view, 
a  +  P  =  7T/4  has  been  shown  to  produce  the  lowest  BER  for  given  values  of  SNR  when 
£  =  0°  [Ref.  8] 
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Figure  5.1       PE  Versus  SNR  for  16-PSK  DBDR  (a  =  11.25°,  e  =  0°). 
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LEGEND 
□    PHS.  ERR  .=  0   PEG. 
o     PHS.  ERR ,=  5   PEG. 
a     PHS.  ERR .=11.25  PEG. 
+     PHS.  ERR/ 
x'"PHSVERRT: 
o    PHS.  ERR.: 

v"phs."err:: 


£2.5  PEG. 
:"3 3.7*5  PEG. 
:  45  PEG. 
="90  PEG." 


a     PHS.  ERR  =135  PEG. 
*    PHS.  ERR  =180  PEG. 


10 


15 


SNR  (DB) 


20 


Figure  5.2    P£  Versus  SNR  for  16-PSK  DBDR  for  Various  Values  of  c  (ct=  11.25°). 
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B.        16-DPSK    DBDR    PERFORMANCE    ANALYSIS    IN    THE    PRESENCE    OF 
RECEIVER  LOCAL  OSCILLATOR  PHASE  ERROR 

Similar  analyses  on  DBDR  performance  have  been  carried  out  for  16-DPSK  (see 

Chapter  IV.B).  so  that  the  notation  introduced  previously  will  be  used  here.   That  is, 

Ck  =  Pr  {  k  correct  }     ,      Ek  =  Pr  {  k  in  error  }  (5.34) 

where  now  k=  1  for  LB,  k=  2  for  MBL,  k=  3  for  MBR  and  k=  4  for  RB. 

In  order  to  evaluate  Pr  {  LB  correct  }  with  differential  decoding,  we  have  to 
evaluate 

Pr{  LB  correct }  =  1/2  [  Pr{  LB  correct  /'  LB  =  0  }  + 

+  {  Pr{  LB  correct  /  LB=  1  }  ]  (5.35) 

where 

Pr{  LB  correct  /  LB  =  0  }  =  (5.36) 

1       1       1 
=  Z    X    £     Pr{  LB  correct/  LB  =  0,  MBL  =  i,  MBR  =  j,  RB  =  k}   x 
i  =  0    j  =  0    k  =  0 

x    Pr{  MBL=i,  MBR  =  j,  RB  =  k} 

where  for  i,  j,  k=  0,  1 

Pr{  MBL  =  i,  MBR  =  j,  RB  =  k}=  1/8  (5.37) 

In  Chapter  IV,  it  has  been  explained  that  with  DPSK  modulation  in  conjuction 
with  bit  direct  detection  methods,  certain  bit  errors  in  the  past  and  presently  recovered 
symbol  states  do  not  result  in  errors  in  the  recovered  data  states.  Such  bit  errors  in 
successive  symbol  states  must  be  considered,  because,  even  though  symbol  errors  may 
occur,  the  delivered  data  may  still  be  correct.  Therefore,  these  errors  have  to  be 
accounted  for  and  considered  in  the  computation  of  the  probability  of  a  bit  being 
correct.  Figure  5.3  shows  one  example  of  the  situation  being  discussed.  That  is,  the 
data  state  (0,0,0,0)  is  assumed  to  occur  and  the  symbol  state  (0,0,0,0)  is  assumed  to 
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CORRECT  RECEPTION  COMBINATIONS  IN  PRESENCE 

OF  ERROR 

nMn»nnnn§nnMMnmmm»»imimM4muum9UumimnimnmnM» 

INPUT           PREVIOUS            PRESENT              PREVIOUS 

PRESENT 

OUTPUT 

TRANSMISSION  TRANSMISSION       RECEPTION 

RECEPTION 

OOOO                0000              0000                0000 

0   0   0   0 

0   0   0   0 

0  0   0   1 

0  0  0   1 

0  0  0   0 

0  0    10 

0   0   10 

0   0   0   0 

0  0    11 

0   0   11 

0   0  0   0 

0   10  0 

0   10   0 

0  0   0  0 

0    10   1 

0    10    1 

0  0   0  0 

0    110 

0   110 

0   0   0   0 

0   111 

0   111 

0  0   0   0 

10  0  0 

10   0   0 

0   0   0  0 

10   0    1 

1  0.  0    1 

0  0   0   0 

10    10 

10   10 

0  0   0   0 

10    11 

10    11 

0   0   0   0 

110   0 

110   0 

0  0   0   0 

110    1 

110    1 

0   0   0   0 

1110 

1110 

0   0   0   0 

1111 

1111 

0   0   0   0 

0   0   0   0 

10   0   0 

0   0    10 

0   0   0    1 

-    1   1   0 

0   0    10 

0   0    10 

110   0 

0   0    10 

0   0    11 

0    111 

0   0    10 

0    10   0 

10    10 

0   0    10 

0   10   1 

0    10   0 

0   0    10 

0    110 

0   0    10 

0   0    10 

0    111 

0   0   0   0 

0   0    10 

10   0   0 

1111 

0   0    10 

10   0    1 

0   0   0    1 

0   0    10 

10    10 

1110 

0   0    10 

10    11 

0    10    1 

0   0    10 

110   0 

110    1 

0   0    10 

110    1 

0   0    11 

0   0    10 

1110 

10   0    1 

0   0    10 

1111 

10    11 

0   0    10 

0   0   0   0 

110    1 

0    10   0 

0   0   0    1 

10    10 

0    10   0 

0   0    10 

10  0    1 

0    10   0 

0   0    11 

0   0    10 

0    10   0 

0    10  0 

1111 

0    10   0 

0    10    1 

10   0   0 

0    10   0 

0   110 

1110 

0    10   0 

0    111 

110   0 

0    10   0 

10   0   0 

0   0    11 

0    10   0 

10   0    1 

0    10   0 

0    10   0 

10    10 

10    11 

0    10   0 

10    11 

0   0   0   0 

0    10   0 

110   0 

0   0   0    1 

0    10   0 

110    1 

0    110 

0    10   0 

1110 

0    10    1 

0    10   0 

1111 

0    111 

0    10   0 

0   0   0   0 

110   0 

0    110 

0   0   0    1 

0    10   0 

0    110 

0   0    10 

1110 

0    110 

0   0    11 

0    110 

0    110 

0    10   0 

10   0   0 

0    110 

0    10    1 

0   0   0   0 

0    110 

0    110 

10    10 

0    110 

0    111 

0   0    10 

0    110 

10   0   0 

110    1 

0    110 

10   0    1 

0    10    1 

0    110 

10    10 

1111 

0    110 

10    11 

0    111 

0    110 

110   0 

10   0    1 

0    110 

110    1 

0   0   0    1 

0    110 

Figure  5.3     Example  of  Error  Patterns  for  16-DPSK  System  Where  Data  State, 
Previous  Transmission,  and  Present  Transmission  are  0000  with  RB  Correctly  Recovered. 
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PRESENT 
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OUTPUT 


Figure  5.3 


Continued. 
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have  been  previously  transmitted.  This  means  that  according  to  the  encoder  matrix 
shown  in  Chapter  II,  the  symbol  state  (0,0,0,0)  must  now  be  generated  for 
transmission.  Figure  5.3  shows  the  receptions  of  the  symbol  states  which  produce 
recovered  data  states  in  which  the  RB  is  correctly  delivered.  Observe  that  there  are  128 
such  combinations  which  yield  0  for  the  RB.  Therefore  notation  has  to  be  introduced 
to  account  for  error  patterns  of  the  above  type,  namely  Pr{  RB  correct,  (g,h,p,q)  / 
LB  =  m,MBL  =  i,MBR  =  j,RB  =  k  },  where  g,h,p,q=  0,1.  The  meaning  of  this 
expression  is  the  probability  of  the  RB  being  correctly  recovered  when  the  data  state  is 
(m,j,i,k)  and  the  delivered  data  state  is  (g,h,p,q)  where  q  =  0  in  this  case  (otherwise  the 
RB  could  not  be  correct  due  to  the  transmission  of  the  RB  =  0).   Therefore 

Pr{  RB  correct  /  LB  =  0,  MBL  =  0  ^BR  =  0,  RB  =  0  }  = 
1       1       1 
=  X      S      £   Pr(  RB  correct,  (g,h,p,0)  /  LB  =  0,  MBL  =  0,  MBR  =  0,  RB  =  0  }  (5.38) 

g  =  0  h  =  0  p  =  0 

Since  each  element  of  the  symbol  state  is  independently  recovered,  we  see  from  Figure 
5.3  that 

Pr{  RB  correct,(0,0,0,0)  /  LB  =  0,MBL  =  0,MBR  =  0,RB  =  0  }  =  CI2  C22  C32  C42  + 
+  CI2  C22  E32  C42  +  CI2  C22  E32  E42+  CI2  E22  C32  E42  +  CI2  E22  E32  C42  + 
+  CI2  E22  E32  E42+  El2  +  C22  C32  C42+  El2  C22  C32  E42+  El2  C22  E32  C42  + 
+  El2  C22  E32  E42+  El2  E22  C32  C42+  El2  E22  E32  C42+  El2  E22  E32  E42  + 
+  CI2  C22  C32  E42  (5.39) 

However,  from  Equations  5.34,  and  5.30,  it  is  shown  that  Cl=  C4  and  El  =  E4. 
Therefore, 

Pr{  RB  correct,  (0,0,0,0)  /  LB  =  0,  MBL  =  0,  MBR  =  0,  RB  =  0}  = 
=  CI4  C22  C32+  2  CI2  C22  C32   El2  +  C22  C32  El2  E22 
+  2  CI4  C22  E32  +  2  CI2   C22  El2  E32+  CI4  C32  E22  +  El4  E22  E32 
+  2  CI2  C32  El2  E22  +  C22  El2   E22  E32  +  CI2  El2  E22  E32  + 
+  C32  El4  E22  +  CI2  C32  El4  (5.40) 
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A  similar  expression  must  be  obtained  for  other  combinations  of  values  of  g,h  and  p  as 
well  as  the  condition  that  RB  is  1.  This  produces  a  total  of  256  possibilities  of  received 
symbol  states  that  result  in  a  correct  RB.  Computer  generated  outputs  were  used  to 
identify  the  equivalent  error  patterns  for  other  cases  involving  LB=1,  MBL=  1  and 
MBR=  1.  These  results  are  shown  in  Appendix  C.  Similar  analyses  must  be  carried 
out  in  order  to  determine  the  unconditional  probabilities  Pr(LB  correct},  Pr{MBL 
correct}  as  well  as  Pr{MBR  correct}.  Such  analyses  demonstrate  due  to  inherent 
symmetries  that 

Pr{LB  correct}  =  Pr{MBL  correct}  =  Pr{MBR  correct}  =  Pr{RB  correct}  =  PE  (5.41) 

so  that  using  Equation  5.38  along  with  similar  forms  obtained  for  k=  1  and  1=  1  yield 

PE  =  1  -  1/2  {  2  CI4  C22  C32  +  CI4  C22  C3  E3  +  3  CI4  C22  C3  E2  + 
+  2  CI4  C2  C3  E2  E3  4-  6  CI3  C22  C32  El  +  4  CI3  C22  C3  El  E3  + 
+  12  CI3  C2  C32  C3  El  E2  +  8  CI3  C2  C3  C3  El  E2  +  6  CI3  C22  El  E32  + 
+  6  CI3  C22  C32  El  +  CI4  C3  E22  E3  +  4  CI3  C3  El  E22  E3  + 
+  2  CI4  El2  E32  +  CI4  C22  E32  +  3  CI4  C2  E2  E32  +  2  C22  C32  El4  + 
+  CI4  C32  E22  +  10  CI2  C22  C32  El2  +  10  CI2  C22  El2  E32  + 
+  10  CI2  C32  El2  E22  +  8  CI2  El2  E22  E32  +  18  CI2  C2  C32  El2  E22  + 
+  18  CI2  C2  El2  E2  E32  +  6  CI3  El  E22  E32  +  12  CI3  C2  El  E2  E32  + 
+  12  CI2  C2  C3  El2  E2  E3  +  6  CI2  C22  C3  El2  E3  +  5  CI2  C3  El2  E22  E3  + 
+  6  CI  C22  C32  El3  +  6  CI  C32  El3  E22  +  6  CI  C22  El3  E22  +  C22  El4  E32 
+  12  CI  C2  C32  El3  E2  +  12  CI  C2  El3  E2  E32  +  6  CI  El3  E22  E32  + 
+  8  CI  C2  C3  El3  E2  E3  +  4  CI  C3  El3  E22  E3  +  4  CI  C22  C3  El3  E3  + 
+  C32  El4  E22+  2  El4  E22  E32  +  C22  C3  El4  E3  +  C3  El4  E22  E3  + 
+  3  C2  El4  E2  E32  +  2  C2  C3  El4  E2  E3  +  2  C2  C32  El4  E2  }  (5.42) 

In  Figure  5.4,  PE  (given  by  Equation  5.42)  is  plotted  for  various  values  of  a.    Also 
Figure  5.5  includes  BER  for  various  values  of  e  when  a=  11.25° 
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C.  NUMERICAL  RESULTS  AND  DISCUSSION  OF  RESULTS  FOR  16-PSK 
AND  16-DPSK  DBDR'S  WITH  RECEIVER  LOCAL  OSCILLATOR  PHASE 
ERROR 

In   this    chapter    16-PSK    and    16-DPSK    DBDRs    have    been   analyzed   when 

operating  in  the  presence  of  a  receiver  local  oscillator  phase  error.    The  performance 

results  which  have  been  obtained  for  the    16-ary  case,  are  similar  to  those  obtain  for 

the  8-ary  case.    Hence,  instead  of  repeating  the  conclusions  similar  to  these  presented 

for  the  8-ary  case,  conclusions  that  are  unique  to  the  16-ary  case  will  be  presented  in 

this  part.  Main  differences  between  8-ary  and  16-ary  cases  which  are  observed  via 

analyses  and  numerical  evaluations  are 

1.  The  optimum  angles  in  minimum  BER  sense,  a  and  P  between  the  symbols 
becomes  11.25°  and  33.75°  respectively  for  16-PSK  case,  as  well  as  for 
16-DSPK  case,  with  a  +  p  =  45° 

2.  The  usage  of  16-PSK  provides  better  performance  than  16-DPSK  for  £^33.5°, 
however  the  actual  BER  values  are  too  high  for  use  in  a  practical  system  when 
e  greater  than  a  few  degrees. 

3.  The  16-DPSK  system  performs  better  than  the  16-PSK  system  for  £>33.5°. 
However  again  the  actual  values  of  PE  under  these  circumstances  is  too  high 
for  a  practical  system. 

4.  The  16-DPSK  system  is  insensitive  to  phase  errors  of  K  radians. 
Performance  comparisons  between  the  16-PSK  and  16-DPSK  DBDRs  are  plotted 

on  Figure  5.6,  and  BERs  are  plotted  on  Figures  5.7  and  5.8  respectively  for  various 
values  of  a  with  £  =  0° 
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Figure  5.6    PE  Versus  SNR  Comparison  Plot  for  16-PSK  and  16-DPSK  DBDR 

(a=  11.25°,  e=  0°). 
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VI.  CONCLUSIONS 

In  this  research  the  use  of  DBDRs  had  been  analyzed  for  M-PSK  and  MDPSK 
communication  systems.  The  main  assumptions  in  the  analyses  are  a  Gray  coding 
mapping  for  the  symbols  and  an  Additive  White  Gaussian  Noise  channel  model.  It  is 
further  assumed  that  the  receiver  has  a  known  phase  error  in  its  local  oscillator  and  its 
effect  on  receiver  performance  is  analyzed  for  the  the  special  cases  in  which  M  =  8  and 
M=16. 

In  addition  to  the  analytical  results,  graphical  results  are  presented  in  the  form  of 
plots  of  variations  in  PE  for  various  values  of  phase  angle  a,  SNR,  and  local  oscillator 
phase  error  £. 

Optimum  conventional  M-PSK  and  M-DPSK  receivers  are  designed  to  detect 
symbols  while  DBDRs  are  designed  to  detect  bits.  However,  local  oscillator  phase 
errors  degrade  the  performance  of  M-PSK  and  M-DPSK  conventional  receivers.  Part 
of  the  work  undertaken  was  intended  to  investigate  the  effect  of  such  phase  errors  on 
DBDR's.  The  variations  of  PE  versus  £  are  plotted  in  Figure  6.1  for  8-PSK  and 
8-DPSK,  and  in  Figure  6.2  for  16-PSK  and  16-DPSK.  If  £  (receiver  local  oscillator 
phase  error)  can  be  constrained  to  a  small  interval,  namely  0°<£<30°,  the  M-PSK 
system  provides  better  performance  than  the  M-DPSK  system.  For  values  of  £  outside 
this  range,  both  systems  have  high  error  rates  that  make  the  receivers  unsuitable  for 
practical  applications.  However,  for  complete  phase  reversal,  i.e.,  £=  180°,  the  M- 
DPSK  system  suffers  no  performance  penalty  whereas  the  M-PSK  system  can  not 
properly  operate. 

The  performance  degradation  due  to  a  phase  error  has  been  analyzed  also  from 
a  different  point  of  view.  That  is,  when  £  *  0°,  there  is  an  increase  in  PE  for  both  the 
M-PSK  and  M-DPSK  receivers.   The  ratio  of  the  phase  error  is  given  by 

PER  =  PE0/  PEj  (6.1) 

where 

PEQ=  Probability  of  bit  error  @  £  =  0° 

?E{  =  Probability  of  bit  error  @  £  *  0° 
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PER  has  been  analyzed  under  various  conditions  for  M-PSK  and  M-DPSK.  The 
results  have  been  plotted  in  Figures  6.3and  6.5  for  M  =  8,  and  in  Figures  6.6-6.8  for 
M=16.  as  a  function  of  SXR  for  different  values  of  e.  The  plots  demonstrate 
significant  differences  in  the  performance  degradation  of  the  M  =  8  versus  M=  16  case. 
However,  for  high  SXR  values,  these  differences  tend  to  vanish. 
The  analysis  results  lead  to  the  following  conclusions 

1.  By  using  DBDRs,  the  error  probability  of  each  bit  in  the  symbol  state  can  be 
calculated  independently  and  more  importantly  these  error  probabilities  can  be 
varied  by  changing  the  angle  between  symbols. 

2.  For  optimum  BER,  a  should  be  set  to  22.5°  for  8-PSK  and  8-DPSK  systems, 
while  16-PSK  and  16-DPSK  systems  require  a=  11.25°  and  P  +  a  =  45°. 

3.  The  solution  for  the  phase  ambiguity  problem  cannot  be  achieved  by  using  M- 
DPSK  DBDR's,  unless  the  phase  ambiguity  is  180°.  Therefore,  if  the  perfect 
coherency  conditions  cannot  be  achieved,  the  M-PSK  system  should  be  used 
whenever  DBDR's  are  used. 

4.  When  e=  0°,  the  DPSK  system  exhibits  good  performance.  However,  if  the 
£=0°  condition  can  be  guaranteed  there  is  no  need  to  use  the  differential 
encoding  scheme. 
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Figure  6.1     PE  Versus  £  for  8-PSK  and  8-DPSK  DBDR  (S\'R=    15  UB  and  a=  22.5°). 
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Figure  6.2    PE  Versus  £  for  16-PSK  and  16-DPSK  DBDR  (SNR=    15dBanda=  11.25°). 
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Figure  6.3     PER  (c  =  0o/c=5°)  Versus  SNR  for  S-PSK 
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Figure  6.5     PER  (s  =  0°/c  =  45°)  Versus  SNR  for  S-PSK 
and  8-DPSK  DBDR's  (a=  22.5°). 
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APPENDIX  A 
DBDR  INPUT/OUTPUT  GENERATION  PROGRAM  (BELIZ) 

1.        USER  GUIDE 

This  program  generates  all  the  possible  input  /  output  combinations  in  an  M- 
DPSK  communication  system.  The  program  runs  interactively  and  sends  outputs  to 
the  files 

1.  FILE  FT08F001 

2.  FILE  RCDATA 

3.  FILE  XMDATA 

4.  FILE  R 

5.  FILE  ERR. 

These  Files  are  used  for  information  transfer  between  subroutines  but  are  not  required 
for  program  execution.  It  is  possible  to  generate  input  /  output  combinations  for 
8-DPSK  and  16-DPSK  by  using  this  program,  while  for  32-DPSK  or  64-DPSK 
systems,  dimensions  and  format  statements  must  be  changed  accordingly,  and  a  larger 
character  variable  set  in  the  subroutine  XMTR  must  be  created.  The  program  accepts 
only  the  integer  equivalents  of  your  data  or  symbol  states.  An  example  run  of  the 
program  is  shown  below 

run  beliz 
EXECUTION  BEGINS... 

PLEASE  SELECT  ONE  OF  THEM 

(1)  IF  YOU  HAVE  TRANSMITTER 

(2)  IF  YOU  HAVE  RECEIVER 


1 
PLEASE  SELECT  ONE  OF  THEM 
:f  YOUR  SYSTEM  8-DPSK 
IF  YOUR  SYSTEM  16-DPSK 


(i)£ 


1 

ENTER  THE  INTEGER  EQUIVALENTS  OF  SYMBOL  STATES 
IN  THE  ORDER  OF  COUNTER  CLOCK  WISE 
ENTER  THE  STATE  #   1 


ENTER  THE  STATE  #  2 
ENTER  THE  STATE  #  3 
ENTER  THE  STATE  #   4 
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3 
2 

0 
4 


ENTER  THE  STATE  #  5 

ENTER  THE  STATE  #  6 

ENTER  THE  STATE  #  7 

ENTER  THE  STATE  #  8 


6 

ENTER  THE  INTEGER  EQUIVALENTS  OF  THE 
PHASE  ANGLES  OF  YOUR  CODING  SYSTEM 
1*ALPHA=  ? 


2*ALPHA=  ? 

1 

3*ALPHA=  ? 

3 

4*ALPHA=  ? 

2 

5*ALPHA=  ? 

6 

6*APLHA=  ? 

7 

7*ALPHA=  ? 

5 

8*ALPHA=  ? 

4 

PLEASE  ENTER  THE  SELECTION  CORRECTLY 

ENTER  1  IF  YOU  NEED  BINARY  FORM  OF  SYMBOL  AND 
DATA  STATES  INCLUDING  SYMBOL  OUTPUTS 

ENTER  2  IF  YOU  NEED  DECIMAL  FORM  OF  SYMBOL  AND 
DATA  STATES  INCLUDING  SYMBOL  OUTPUTS 

ENTER  3  IF  YOU  WANT  TO  TURN  BACK  MAIN  PROGRAM. 

1 

PLEASE  ENTER  THE  SELECTION  CORRECTLY 

ENTER  1  IF  YOU  NEED  BINARY  FORM  OF  SYMBOL  AND 
DATA  STATES  INCLUDING  SYMBOL  OUTPUTS 

ENTER  2  IF  YOU  NEED  DECIMAL  FORM  OF  SYMBOL  AND 
DATA  STATES  INCLUDING  SYMBOL  OUTPUTS 

ENTER  3  IF  YOU  WANT  TO  TURN  BACK  MAIN  PROGRAM. 

3 

ANALYSES  RESULTS  WAS  SENT  TO  THE  FILE  FT08F001 
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SOURCE  CODE 


c 

C  THIS  PROGRAM   SIMULATES  A  M-DPSK  COMMUNICATION  SYSTEM 

C  PROGRAM  RUNS  TOTATLY  INTERACTIVE  THEN  SENDS  TO  OUTPUTS  FILES; 

C  FILE  FT08F001,FILE  RCDATA,FILE  XMDATA,FILE  R,FILE  ERR. 

C  THESE  FILES  ARE  USED  FOR  INFORMATION  TRANSFER  BETWEEN 

C  SUBROUTINES  BUT  ARE  NOT  REQUIRED  FOR  PROGRAM  EXECUTION. 

C  YOU  CAN  SIMULATE  8-DPSK  AND  16-DPSK  BY  USING  THIS  PROGRAM 

C  FOR  32-DPSK  OR  64-DPSK  SYSTEMS  YOU  HAVE  TO  CHANGE  DIMENSIONS 

C  AND  FORMAT  STATEMENTS  ACCORDINGLY, ALSO  YOU  HAVE  TO  EXPEND 

C  CHARACTER  VARIABLE  SET  IN  THE  SUBROUTINE  XMTR. 

C 

C  OPERATION  SPECIFICATIONS ; 

C 

C  THE  SYMBOL  STATES  REQUESTEd  DEPEND  ON  YOUR  SYSTEM  AND  CAN 

C  BE  ANY  SEQUENCE.   MAIN  THING  IS;  YOU  HAVE  TO  SELECT  YOUR  REFERENCE 

C  THEN  ENTER  THEM  BY  ONE  STEP  INCREMENTS , ALSO  SAME  SPECIFICATION 

C  EXISTS  FOR  YOUR  ANGLES .PROGRAM  ONLY  EXCEPTS  THE  INTEGER  EQUIVALENTS 

C  OF  YOUR  STATE  OR  ANGLE  CODES. 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 
c 

C      VARIABLE  DECLERATION 
C 

INTEGER  N,L, ORDER, TM(64, 64) ,PS(64) ,TRDATA(256) ,M,KK,K,THETA(64) 

INTEGER  MAT(64) , KL ,RMAT(64, 64) ,TR(64) ,TH(64 , 64) ,RCVMAT(64, 64) ,NN 

INTEGER  QQU,QU,SQU,LD,IS.NT(1,8) .KM(1,8),KN(1,8) 

INTEGER  RCV(64,64) ,T(256) ,TD(256) ,NEXTST(256) ,PRESTA(256) ,FL 

INTEGER  A ,  B  ,  Y ,  KD  ,  KY  ( 1 ,  8 )  ,  KH  ( 1 ,  8 )  ,  K  J  ( 1 ,  8 )  ,  KB  ( 1 ,  8 ) 

CHARACTERS  DATA (64, 64) 

CHARACTER*16  SET 
C 
C      ***     CALLING  SUBROUTINES      *** 

CALL  INFO  (ORDER, PS, THETA,NN,N,FL) 

IF  (FL.EQ.l)THEN 

CALL  XMTR  (ORDER, FL) 

CALL  XMSTA  ( THETA,TM, PS , ORDER) 

CALL  XMDAT ( TM , ORDER , THETA , NEXTST , PS , T , TD ) 

CALL  CONVER  ( N, NEXTST ,NN, ORDER, FL) 

WRITE(6,1) 

1  FORMAT(T9,'  ANALYSES  RESULTS  WAS  SENT  TO  THE  FILE  FT08F001 ■ ) 

ELSE 

CALL  RCVR( ORDER) 

CALL  RCVSTA  (PS , ORDER, THETA, RCVMAT) 

CALL  RCVDAT ( TR , PS , ORDER , NEXTST , RCVMAT ) 

CALL  CONVER  ( N, NEXTST, NN, ORDER, FL) 

CALL  ERROR ( FL , THETA , RCVMAT , PS , ORDER , TR , TM , T , NEXTST , TD ) 

WRITE(6,2) 

2  FORMAT(T9,'  ANALYSES  RESULTS  WAS  SENT  TO  THE  FILE  FT08F001',/, 
*T7,' ERROR  ANALYSES  RESULTS  WAS  SENT  TO  THE  FILE  ERR1) 

END  IF 

STOP 

END 
C 
C 

SUBROUTINE  INFO ( ORDER , PS , THETA , NN , N , FL ) 
C 

C    THIS  SUBROUTINE  PROVIDES  INFORMATIONS  ABOUT  OUR  SYSTEM 
C 

INTEGER  N.L, ORDER, TM(64, 64) ,PS(64) ,THETA(64) , FL 

CHARACTERS  DATA (64, 64) 

WRITE  (6,*)  '  PLEASE  SELECT  ONE  OF  THEM' 

WRITE  (6,*)  '(1)  IF  YOU  HAVE  TRANSMITTER' 

WRITE  (6 ,*)  '(2)  IF  YOU  HAVE  RECEIVER1 

READ  (6,*)  FL 

WRITE  (6,*)  ■  PLEASE  SELECT  ONE  OF  THEM' 

WRITE  (6,*)  '(1)  IF  YOUR  SYSTEM  8-DPSK' 
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WRITE  (6 .*)  '(2)  IF  YOUR  SYSTEM  16-DPSK' 

READ  (6,*)  L 

IF(L.EO.l)  THEN 

ORDER=8" 

NN=63 

N=6 

END  IF 

IF(L.EQ.2)  THEN 

ORDER=16 

N=8 

NN=255 

ENDIF 

WRITE  (6,*)  'ENTER  THE  INTEGER  EQUIVALENTS  OF  SYMBOL  STATES' 

WRITE  (6,*)  'IN  THE  ORDER  OF  COUNTER  CLOCK  WISE' 

DO  11  K  =1. ORDER 

WRITE (6 ,20)   K 
20   FORMAT (1 OX, 'ENTER  THE  STATE  #',I3) 

READ(6,*)  PS(K) 
11   CONTINUE 

WRITE  (6,*)  'ENTER  THE  INTEGER  EQUIVALENTS  OF  THE  ' 

WRITE  (6,*)  'PHASE  ANGLES  OF  YOUR  CODING  SYSTEM  ' 

DO  33  L=l, ORDER 

WRITE  (6,10)  L 
10    FORMAT  (10X,I3, ' *ALPHA=  ?') 

READ(6,*)  THETA(L) 
33    CONTINUE 

RETURN 

END 
C 

SUBROUTINE  XMTR ( ORDER, FL) 
C 
C 

C        THIS  SUBROUTINE  PROVIDES  GENERAL  TRANSMITER  MATRIX  FOR 
C     M-DPSK  SYSTEM. 

C        ITS  OUTPUT  CAN  BE  SEEN  ON  EITHER  FILE  FT08F001  OR  FILE 
C     XMDATA. 

INTEGER  N.L, ORDER, FL 

CHARACTER*16  SET 

SET(1:16)=  ' ABCDEFGHIJKLMNOP ' 

OPEN ( UNI T= 1 1,F I LE=' XMDATA' , STATUS= ' OLD ' ) 

IF  (FL.EQ.l)  THEN 

WRITE  (8,81)0RDER 
81    FORMAT (///,T1 5, 'TRANSMITTER'S  ENCODER  MATRIX  ANALYSIS 
*  FOR' ,13, '-DPSK' 
*,/,  60 ('-'),//, T5 , ' GENERAL  TRANSMITTER  MATRIX ' , / , T2 , 31 ( ' - ■ ) ) 

ELSE 

WRITE  (8,88) 
88   FORMAT (///,T1 5, 'TRANSMITTER  INFORMATION  FOR  THE  GIVEN  RECEIVER' 
*,/,72('-'),//,T5, 'GENERAL  TRANSMITTER  MATRIX' ,/ ,T2, 31 ('-') ) 

ENDIF 

DO  3  N=l, ORDER 

IF  (N.EQ.l)  THEN 

WRITE  (§,*)  SET(N:ORDER) 

WRITE  (11,*)  SET(N:ORDER) 

ELSE 

WRITE  (8,*)  SET(NtORDER) ,SET(1:N-1) 

WRITE  (11,*)  SET(N:ORDER),SET(l:N-l) 

ENDIF 
3    CONTINUE 

CLOSE  (11) 

RETURN 

END 
C 

SUBROUTINE  XMSTA ( THETA , TM , PS , ORDER ) 
C 

C     THIS  SUBROUTINE  CONVERTS  GENERAL  XMTR  MATRIX  IN  TO  OUR  SPECIAL 
C    XMTR  SYMBOL  MATRIX. 

C     GENERAL  XMTR  MATRIX  IS  READ  FROM  FILE  XMDATA. 
C 

INTEGER  TM ( 64 , 64 ) , PS ( 64 ) , THETA ( 64 ) , ORDER 
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CHARACTER*!  DATA (64, 64) 

OPEN ( UNIT=1 1 , FILE= ' XMDATA ' , STATUS  = ' OLD ' ) 

IF  ( ORDER. EQ. 8) THEN 

WRITE(8,99)(THETA(I) ,1=1, ORDER) 
99    FORMAT (///,T5, 'TRANSMITTER  MATRIX  ',/, 31 ('-'),// , 
*' SYMBOL  STATE  (COLUMNS)  /  DATA  STATE  (ROWS)' 
*,//,4X,8(I3K/,35('-')K 

READ  (11,10)  ((DATA(L,N),N=1,8),L=1,8) 

10  FORMAT  (1X,8A1) 

ELSE 

WRITE (8, 43 )(THETA(I), 1=1, ORDER) 
43    FORMAT (///,T5, 'TRANSMITTER  MATRIX  ',/, 31 ('-'),// , 
*' SYMBOL  STATE  (COLUMNS)  /  DATA  STATE  (ROWS)' 
*,//,4X, 16(13),/, 70('-'  ) 

READ  (11,20)  ((DATA(L,N),N=1,16),L=1,16) 
20    FORMAT  (1X,16A1) 

END  IF 

DO  3  K=l, ORDER 

DO  4  N=l, ORDER 

DO  6  M=  0, ORDER 

IF(ICHAR(DATA(K,N)).EQ.ICHAR(DATA(l,ORDER-M)))THEN 

TM(K,N)=PS(ORDER-M) 

ENDIF 

6  CONTINUE 

4  CONTINUE 

IF  ( ORDER. EQ. 8) THEN 
WRITE(11,11)  (TM(K,N),N=1. ORDER) 
WRITE(3,83)  PS(K).(TM(K,N)/N=1, ORDER) 
83   FORMAT (13,'  ',8(13)) 

11  FORMAT  (812) 

ELSE 
77    WRITE(11,12)  (TM(K,N),N=1, ORDER) 

WRITE(8,93)  PS(K) , (TM(K,N) ,N=1, ORDER) 
93    FORMAT(I3, ' I ' ,16(13)) 

12  FORMAT  (1613) 

ENDIF 
3   CONTINUE 

RETURN 

END 
C 

SUBROUTINE  XMDAT ( TM , ORDER , THETA , NEXTST , PS , T , TD ) 
C 

C         THIS  SUBROUTINE  FINDS  THE  NEXT  SYMBOL  STATE/S  BY  DETERMINING 
C     XMTR  MATRIX  CONSTRUCTED  BY  SUB.  XMSTA 
C 
C 

INTEGER  TM(64,64) ,PS(64) ,TRDATA(256) ,KK.K,THETA(64) ,NN 

INTEGER  ORDER, T(256) ,TD(256) ,NEXTST(256) ,PRESTA(256) 

DO  5  1=1, ORDER 

T(PS(I))=I 

TD(THETA(I))=I 

5  CONTINUE 

KK=0 

DO  7  K=0, ORDER- 1 
DO  8  M=0, ORDER- 1 
KK=KK+1 

MEXTST(KK)=TM(T(K),TD(M)) 
8    CONTINUE 

7  CONTINUE 

RETURN 

END 
C 

SUBROUTINE  RCVR  (ORDER) 
C 

C  THIS  SUBROUTINE  PREPARES  GENERAL  RECEIVER  MATRIX  AND 

C     SENDS  OUTPUTS  TO  THE  "FILE  RCDATA" . 
C 

INTEGER   ORDER , MAT ( 64 ) , KL , RMAT (64,64) 

OPEN ( UNIT=2 . FILE= ' RCDATA ■ , STATUS= ' OLD ' ) 

WRITE  (8,31)ORDER 
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81    F0RMAT(///,T15, 'RECEIVER'S  DECODER  MATRIX  ANALYSIS 
*  FOR' ,13, ' -DPSK' ,/, 

*60 ('-'), //,T5, 'GENERAL  RECEIVER  MATRIX' ,/ ,T2, 31 ('-') ) 
DO  3  1=0, ORDER- 1 

MAT(I+1)=I 
3    CONTINUE 

DO  44  K=l, ORDER 

RMAT(1,K)=MAT(K) 
44    CONTINUE 

DO  55  L=2, ORDER 

DO  66  M=1,L-1 

RMAT(L,M)=MAT(0RDER-L+M+1) 
66    CONTINUE 

DO  77  N=L, ORDER 

RMAT(L,N)=MAT(N-L+1) 
77    CONTINUE 
55   CONTINUE 

DO  88  NK=1, ORDER 

IF  (ORDER. EQ. 3) THEN 

WRITE  (2,10)(RMAT(NK,KL) ,KL=1, ORDER) 

WRITE  (8,15) (RMAT(NK,KL) ,KL=1, ORDER) 

10  F0RMAT(8(I2)) 

15  F0RMAT(8X,3(I2)) 

ELSE 

WRITE  (2,11)(RMAT(NK,KL) ,KL=1, ORDER) 

WRITE  ( 8 , 16 ) ( RMAT (NK , KL ) , KL=1 , ORDER) 

11  FORMAT (16 (12)) 

16  FORMAT (8X, 16(12)) 

END  IF 
88    CONTINUE 

CL0SE(2) 

RETURN 

END 
C 

SUBROUTINE   RCVSTA  (PS , ORDER, THETA , RCVMAT ) 
C 

C        THIS  SUBROUTINE  PREPARES  SPECIFIC  RECEIVER  MATRIX  UNDER 
C     DEFINED  STATES. IT  READS  GENERAL  RECEIVER  MATRIX  FORM  FILE 
C     RCDATA  WHICH  HAD  BEEN  CREATED  BY  THE  SUBROUTINE  RCVDAT. 
C       OUTPUT  ARE  SENT  TO  THE  FILE  R. 
C 

INTEGER  TH ( 64 , 64 ) , ORDER , PS ( 64 ) , THETA ( 64 ) , RCVMAT (64,64) 

OPEN (UNIT=2,FILE=' RCDATA' , STATUS= ' OLD ■ ) 

OPEN(UNIT=20,FILE='R' , STATUS= ' OLD ' ) 

IF  ( ORDER. EQ. 8) THEN 

WRITE(8, 99) (PS(L),L=1, ORDER), (THETA(I), 1=1. ORDER) 
99    FORMAT (///,T5, 'RECEIVER  MATRIX  IS ',/, 31 ('-'),// . 

*' SYMBOL  STATE  (COLUMNS)  /  DATA  STATE  (ROWS) ',// ,4X, 8(13) ) 

ELSE 

WRITE(8, 199) (PS(L),L=1, ORDER) 
199    FORMAT (///,T5, 'RECEIVER  MATRIX  IS ',/, 31 ('-'),// , 
*' SYMBOL  STATE  (COLUMNS)  /  DATA  STATE  (ROWS)' 
*,//,4X, 16(13)) 

END  IF 

DO  1  1=1, ORDER 

IF  ( ORDER. EQ. 8) THEN 

READ(2,11)  (TH(I,N),N=1, ORDER) 

11  FORMAT (812) 

ELSE 

READ(2,12)  (TH(I,N),N=1, ORDER) 

12  FORMAT (1612) 

END  IF 

DO  2  K=l, ORDER 
RCVMAT(I,K)=THETA(TH(I,K)+1) 
2    CONTINUE 

IF  (ORDER. EQ.8)THEN 

WRITE(20,13)  (RCVMAT(I,N) ,N=1, ORDER) 

13  FORMAT (812) 

WRITE(8,83)  PS(I).(RCVMAT(I,N),N=1, ORDER) 
83    FORMAT ( 13, ' I ' ,8(13)) 
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ELSE 

WRITE(8,93)  PS(I),(RCVMAT(I,N),N=1, ORDER) 

93  FORMAT (13, ' f  \  16(13) ) 

WRITE (20,14)  (RCVMAT ( I , N) , N=l , ORDER) 
14    FORMAT(16I2) 
END  IF 

1  CONTINUE 

CLOSE(20) 

RETURN 

END 
C 

SUBROUTINE  RCVDAT (TR , PS , ORDER , NEXTST , RCVMAT) 
C 

C         THIS  SUBROUTINE  READS  SPECIFIC  RECEIVER  MATRIX  WHICH  HAD 
C       BEEN  CREATED  BY  SUBROUTINE  RCVSTA,THEN  SEARCHS  FOR  THE  RCVR 
C       OUTPUT. 
C 
C 

INTEGER  RCVMAT ( 64 , 64 ) , PS ( 64 ) , KK , ORDER 

INTEGER  TR(64),NEXTST(256) 

OPEN(UNIT=20,FILE='R' , STATUS= ' OLD ' ) 

DO  23  LL=1, ORDER 

IF  ( ORDER. EO. 8) THEN 

READ  (20 , 24 J (RCVMAT (LL , KM) , KM=1 , ORDER) 
24   FORMAT(8(I2)) 

ELSE 

READ  (20,94)(RCVMAT(LL,KM),KM=1,ORDER) 

94  FORMAT (16 (12)) 

END  IF 
23   CONTINUE 

DO  2  1=1, ORDER 

M=PS(I) 

TR(M)=I 

2  CONTINUE 

KK=0 

DO  3  M=0, ORDER- 1 
DO  4  K=0,ORDER-1 
KK=KK+1 

NEXTST(KK)=RCVMAT(TR(M),TR(K)) 
4   CONTINUE 

3  CONTINUE 

RETURN 
END 


ENTER  1  IF  YOU  NEED  BINARY  FORM  OF  SYMBOL  AND 
DATA  STATES  INCLUDING  SYMBOL  OUTPUTS' 

ENTER  2  IF  YOU  NEED  DECIMAL  FORM  OF  SYMBOL  AND 
DATA  STATES  INCLUDING  SYMBOL  OUTPUTS' 

ENTER  3  IF  YOU  WANT  TO  TURN  BACK  MAIN  PROGRAM. ' 


SUBROUTINE  CONVER (N, NEXTST, NN, ORDER, FL) 

INTEGER  N,S,QU,SU,LLL,LLK(50) ,MMM,NN, FL ,MM,M,K(9 , 1024) 
INTEGER  MW,NT(256,4) ,KZ, ORDER, NEXTST(256) ,ST,QQU,KD 

300  WRITE  (6,*)   'PLEASE  ENTER  THE  SELECTION  CORRECTLY' 

WRITE  (6,** 

WRITE  (6,* 

WRITE  (6,* 

WRITE  (6,* 

WRITE  (6  * 

READ  (6,*)  S 

IF  (S.EQ.l)  GOTO  100 

IF  (S.EQ.2)  GOTO  200 

IF  (S.EQ.3)  GOTO  400 

IF  (S.NE.1.0R.S.NE.2.0R.S.NE.3)  GOTO  300 
200   IF  (FL.EQ.l)  THEN 

WRITE(8,6) 
6   FORMAT  (//,8X,'DATA  STATE ', 5X, ' SYMBOL  STATE ', 5X, ' SYMBOL ' ) 

ELSE 

WRITE(8,21) 
21   FORMAT  (//,3X, 'DATA  STATE' ,5X, 'SYMBOL  STATE1 ,5X, 'SYMBOL') 

ENDIF 

KK=0 
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DO  7  KZ=0, ORDER- 1 

DO  8  M=0, ORDER- 1 

KK=KK+1 

WRITE(8,9)   KZ,M,  NEXTST(KK) 
9    FORMAT  (10X,I3,15X,I3,15X,I3) 
8   CONTINUE 
7    CONTINUE 

GOTO  300 
100   IF  (FL.EQ.l)  THEN 

WRITE (8, 47) 
47   FORMAT  (//.7X,'DATA  STATE ' , 10X, ' SYMBOL  STATE', 1 OX, 
*' SYMBOL' ,f,12{x-x)) 

ELSE 

WRITE (3 ,77) 
77   FORMAT  (// ,7X, ■ SYMBOL  STATE ', 10X, 'DATA  STATE', 10X, 
*' SYMBOL' ,/, 72 ('-')) 

ENDIF 

IF(0RDER.EQ.8)THEN 

KD=3 

ELSE 

KD=4 

ENDIF 

DO  44  L=1,NN+1 

QOU=NEXTST(L) 

DO  33  SR=KD,1,-1 

NT(L,SR)=M0D(QQU,2) 

QQU=QQU/2 
33   CONTINUE 
44   CONTINUE 

DO  23  MW=0,NN 

QU=MW 

DO  22  LLL=N,1,-1 

K(MW,LLL)=M0D(QU,2) 

QU=QU/2 

22  CONTINUE 

IF  ( ORDER. EQ. 8) THEN 

WRITE (8, 45)  MW,(K(MW,LK),LK=1,N),(NT(MW+1,ST),ST=1,KD) 

45  FORMAT  (13 , 2X,3(2X, 12 ,2X) , 3X,3(2X, 12, 2X) ,3X,3(2X, 12 , 2X) ) 

ELSE 

WRITE(8,55)  MW, (K(MW,LK) ,LK=1 ,N) , (NT(MW+1 ,ST) ,ST=1 ,KD) 
55    FORMAT  (13 ,2X,4(2X, 12, IX) ,3X,4(2X, 12 , IX) ,3X,4(2X, 12, IX) ) 
ENDIF 
WRITE  (8,46) 

46  FORMAT  (72('-')) 

23  CONTINUE 

GOTO  300 
400   RETURN 

END 
C 

SUBROUTINE  ERROR (FL,THETA,RCVMAT, PS , ORDER, TR , TM , T , NEXTST , TD ) 
C 

C      THIS  SUBROUTINE  LOCATES  UNDETECTED  ERROR  COMBINATIONS. 
C 
C 

INTEGER  QQU,QU. SOU, PS(64), ORDER, RCVMAT(64, 64), N,TR(64) ,TM(64,64) 

INTEGER  TD(256),T(256),A,B,Y.KD,KN(1,8),KJ(1,8),KB(1,8),KY(1,8) 

INTEGER  FL,IS,NT(1,8),KM(1,8),NEXTST(256),THETA(64),KH(1,8) 

OPEN(UNIT=l , FILE= ' ERR ' , STATUS= ' OLD ' ) 

OPEN(UNIT=15,FILE='ERDATA' , STATUS= ' OLD ' ) 

WRITE  (6 ,*)  '    IF  YOU  WANT  TO  HAVE  UNDETECTABLE  ERROR' 

WRITE (6,*)  'ANALYSES  ENTER  "1" , OTHERWISE  ENTER  "2" ' 

READ  (6,*)  LD 

IF  ( ORDER. EQ. 8) THEN 

KD=3 

ELSE 

KD=4 

ENDIF 

IF(LD.EQ.1)THEN 

WRITE  (b,*)'    WHAT  KIND  OF  ANALYSES  OUTPUT  DO  YOU  WANT' 


it: 


WRITE  (6,*)'   (1)  IN  BINARY  FORM' 
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WRITE  (6 ,*)'   (2)  IN  DECIMAL  FORM1 

READ  (6,*)  IS 

WRITE (1,13) 
13    FORMAT(////,72('#' ),//,7X, 'CORRECT  RECEPTION  COMBINATIONS  IN  PRESE 
*NCE  OF  ERROR' ,//, 72 ('#' )) 

WRITE (1,12) 
12   FORMAT (//,T5, 'INPUT' ,T17, 'PREVIOUS' ,T30, 'PRESENT' ,T43, 'PREVIOUS' , 
*T55, 'PRESENT' ,T63, 'OUTPUT' ,/,T15, 'TRANSMISSION' ,T28, 'TRANSMISSION' 
* , T42 , ' RECEPTION '  , T54 , ' RECEPTION '  , /  ,  79 ( ' - ' ) ) 

ELSE 

GOTO  111 

END  IF 

FL=2 

CALL  XMTR  (ORDER, FL) 

CALL  XMSTA  ( THETA,TM, PS , ORDER) 

CALL  XMDAT(TM, ORDER, THETA,NEXTST, PS ,T,TD) 

DO  1  I=0,ORDER-l 

DO  2  K=0. ORDER- 1 

A=TM(T(I) ,TD(K)) 

B=RCVMAT(TR(I),TR(A)) 

IF(IS.EQ.1)THEN 

OQU=A 

DO  333  JJ=KD,1,-1 

NT(l,JJ)=MOD(QQU,2) 

QQU=QOU/2 
333   CONTINUE 

QU=K 

DO  222  LM=KD,1,-1 

KM(l,LM)=MOD(0_U,2) 

QU=QU/2 
222   CONTINUE 

SQU=I 

DO  444  LN=KD,1,-1 

KN(l,LN)=MOD(SQU,2) 

SQU=SQU/2 
444   CONTINUE 

SQU=B 

DO  555  LN=KD,1,-1 

KB(l,LN)=MOD(SQU,2) 

SQU=SOU/2 
555   CONTINUE 

IF  ( ORDER. EQ. 8) THEN 

WRITE(1,45)(KM(1,L),L=1,KD),(KN(1,KK),KK=1,KD),(NT(1,IT),IT=1,KD) 
*,(KN(l,KK),KK=l,KD),rNT(l,IT),IT=l,KD),(KB(l/IJ),IJ=l.KD) 
45    FORMAT(T3,3(I3),T15,3(I3),T28,3(I3),T42,3(I3),T54,3(I3),T65,3(I3)) 

WRITE ( 1 5 , 85 ) (KN ( 1 , KK) , KK=1 , KD ) , (NT ( 1 , IT) , IT=1 , KD ) 
*,(KB(l.IJ).li=l,KD) 
85   FORlJIAT(9(I2)) 

ELSE 

WRITE(1,55)(KM(1,L),L=1,KD),(KN(1,KK),KK=1,KD),(NT(1,IT),IT=1,KD) 
*,  (K1«I(1,KK)  ,KK=1,KD)  ,  (NT(1,IT)  ,IT=1,KD)  ,  (KB(1,IJ)  ,IJ=1.KD) 
55   FORMAT (Tl, 4( 12 ), T15, 4(12 ), T28,4( 12 ),T42,4( 12), T54.4( 12), T65, 4(12)) 

WRITE(15,95)(KN(1,KK),KK=1,KD),(NT(1,IT),IT=1,KD) 

*,(KB(1,IJ),IJ=1,KD) 
95   FORMAT (12 (12)) 

END  IF 

ELSE 

WRITE (1,10)  K,I,A,I,A,B 
10   FORMAT(T10,I3,T25,I3,T35,I3,T48,I3,T58,I3,T71,I3) 

END  IF 

DO  3  L=0, ORDER- 1 

DO  4  J=0,ORDER-1 

IF(L.EQ.I.AND.J.EQ.A)THEN 

GOTO  100 

ELSE 

Y=RCVMAT(TR(L) ,TR( J) ) 

IF(Y.EQ.B)THEN 

IF(IS.EQ.1)THEN 

SQU=L 

DO  666  LN=KD,1,-1 
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KH(1,LN)=M0D(SQU,2) 

S0U=SQU/2 
666   CONTINUE 

SOU=J 

DO  777  LN=KD,1,-1 

KJ(l,LN)=MOD(SQU,2) 

SQU=SOU/2 
777   CONTINUE 

SOU=Y 

DO  888  LN=KD,1,-1 

KY(l,LN)=MOD(SQU,2) 

SQU=SQU/2 
888   CONTINUE 

IF(ORDER.EQ.8)THEN 

WRITE(1,48J(KH(1,KK) ,KK=1,KD) , (KJ(1,KK) ,KK=1,KD) 
*,(KY(1,IT),IT=1.KD) 

48  FORMAT  (T42 , 3(13) ,T54, 3(13) ,T65.3(I3) ) 

WRITE ( 1 5 , 9 9 9  J t  ( KH ( 1 , KK ) , KK= 1 , KD ) , ( K J ( 1 , KK ) , KK= 1 , KD ) 
*,(KY(1,IT) ,IT=1,KD) 
999    FORMAT  (9(12)) 

ELSE 

WRITE(1,49)(KH(1/KK)/KK=1,KD),(KJ(1;KK),KK=1/KD) 
*,(KY(1,IT),IT=1.KD) 

49  FORMAT  (T42 ,4(12) ,T54, 4(12) ,T65 ,4(12) ) 

WRITE (15, 949 )(KH(1,KK),KK=1,KD),(KJ(1/KK),KK=1/KD) 
*,(KY(1,IT),IT=1,KD) 
949    FORI^IAT  (12(12)) 

END  IF 

ELSE 

WRITE (1,20)  L,J,Y 
20    FORMAT(T48,I3,T58,I3,T71,I3) 

END  IF 

END  IF 
100   ENDIF 
4    CONTINUE 
3     CONTINUE 

WRITE  (1,46) 
46    FORMAT  (79 ('-■)) 
2    CONTINUE 
1    CONTINUE 
111   CLOSE(l) 

CLOSE (15) 

RETURN 

END 
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APPENDIX  B 
COMPUTER  GENERATED  ENCODER  AND  DECODER  OUTPUTS 

1.       8-DPSK  TRANSiMITTER 

TRANSMITTER'S   ENCODER      MATRIX  TABLE   FOR      8-DPSK 
GENERAL   TRANSMITTER  MATRIX 


ABCDEFGH 
BCDEFGHA 
CDEFGHAB 
DEFGHABC 
EFGHABCD 
FGHA3CDE 
GHABCDEF 
HA3CDEFG 


TRANSMITTER  MATRIX 


SYMBOL   STATE    (COLUMNS)    /    DATA   STATE    (ROWS) 
/  01326754 


75132046 
51320467 
13204675 
32046751 
20467513 
04675132 
46751320 
67513204 


SYMBOL 

STATE 

DATA  STATE 

SYMBOL 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

1 

0 

0 

2 

0 

0 

0 

0 

1 

0 

1 

1 

1 

3 

0 

0 

0 

0 

1 

1 

1 

1 

0 

4 

0 

0 

0 

1 

0 

0 

0 

1 

0 

5 

0 

0 

0 

1 

0 

1 

0 

1 

1 

6 

0 

0 

0 

1 

1 

0 

1 

0 

1 

7 

0 

0 

0 

1 

1 

1 

0 

0 

1 

8 

0 

0 

1 

0 

0 

0 

0 

0 

1 

9 

0 

0 

1 

0 

0 

1 

0 

1 

1 

10 

0 

0 

1 

0 

1 

0 

0 

0 

0 

11 

0 

0 

1 

0 

1 

1 

0 

1 

0 

12 

0 

0 

1 

1 

0 

0 

1 

0 

1 
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13 

0 

0 

1 

1 

0 

1 

1 

1 

1 

14 

0 

0 

1 

1. 

1 

0 

1 

0 

0 

15 

0 

0 

1 

1 

1 

1 

1 

1 

0 

16 

0 

1 

0 

0 

0 

0 

0 

1 

0 

17 

0 

1 

0 

0 

0 

1 

0 

0 

0 

18 

0 

1 

0 

0 

1 

0 

1 

1 

0 

19 

0 

1 

0 

0 

1 

1 

1 

0 

0 

20 

0 

1 

0 

1 

0 

0 

0 

1 

1 

21 

0 

1 

0 

1 

0 

1 

0 

0 

1 

22 

0 

1 

0 

1 

1 

0 

1 

1 

1 

23 

0 

1 

0 

1 

1 

1 

1 

0 

1 

24 

0 

1 

1 

0 

0 

0 

0 

1 

1 

25 

0 

1 

1 

0 

0 

1 

0 

1 

0 

26 

0 

1 

1 

0 

1 

0 

1 

0 

0 

27 

0 

1 

1 

0 

1 

1 

0 

0 

0 

28 

0 

1 

1 

1 

0 

0 

0 

0 

1 

29 

0 

1 

1 

1 

0 

1 

1 

0 

1 

30 

0 

1 

1 

1 

1 

0 

1 

1 

0 

31 

0 

1 

1 

1 

1 

1 

1 

1 

1 

32 

1 

0 

0 

0 

0 

0 

1 

0 

0 

33 

1 

0 

0 

0 

0 

1 

1 

1 

0 

34 

1 

0 

0 

0 

1 

0 

1 

0 

1 

35 

1 

0 

0 

0 

1 
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2.        8-DPSK  RECEIVER 


RECEIVER'S   DECODER  MATRIX  TABLE   FOR      8-DPSK 

GENERAL   RECEIVER  MATRIX 
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12345670 
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TRANSMITTER 

INFORMATION 

FOR  THE 

GIVEN 

RECEIVER 

GENERAL  TRANSMITTER  MATRIX 


ABCDEFGH 
BCDEFGHA 
CDEFGHAB 
DEFGHABC 
EFGHABCD 
FGHABCDE 
GHABCDEF 
HABCDEFG 


TRANSMITTER  MATRIX 


SYMBOL  STATE  (COLUMNS)  /  DATA  STATE  (ROWS) 
01326754 


7  75132046 

5  51320467 

1  13204675 

3  32046751 

2  20467513 
0  04675132 

4  46751320 

6  67513204 


102 


3.        16-DPSK  TRANSMITTER 


TRANSMITTER'S  ENCODER  MATRIX  TABLE  FOR   8-DPSK 


GENERAL  TRANSMITTER  MATRIX 


ABCDEFGHIJKLMNOP 
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16-DPSK  RECEIVER 


RECEIVER'S  DECODER  MATRIX  TABLE  FOR   16-DPSK 


GENERAL  RECEIVER  MATRIX 
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APPENDIX  C 

COMPUTER  GENERATED  OUTPUTS  FOR  DBDR'S  ERROR  PATTERN 

ANALYSES 

The  program  outputs  given  below  are  a  partial  sample  of  the  main  outputs. 

Since  undetectable  error  patterns  are  similar,  in  the  next  sections,  1/8  period  of  the 

pattern  have  been  included.  Actually  526  combinations  of  error  patterns  are  detected 
by  the  program  for  the  8-DPSK  case.  For  the  16-DPSK  case  there  are  4096 
undetectable  error  pattern  combinations. 

1.        8-DPSK  RECEIVER 

CORRECT  RECEPTION  COMBINATIONS  IN  PRESENCE  OF  ERROR 
###################################################################### 

INPUT     PREVIOUS     PRESENT      PREVIOUS  PRESENT  OUTPUT 

TRANSMISSION  TRANSMISSION   RECEPTION  RECEPTION 

000  000      000       000  000  000 

001  001  000 

010  010  000 

011  011  000 

100  100  000 

101  101  000 

110  110  000 

111  111  000 

001  000      100       000  100  001 

001  Oil  001 

010  000  001 

Oil  010  001 

100  110  001 

101  001  001 

110  111  001 

111  101  001 

010     000      111       000  111  010 

001  000  010 

010  110  010 

Oil  100  010 

100  101  010 

101  010  010 

110  001  010 

111  Oil  010 

Oil     000      110       000  110  Oil 

001  010  Oil 

010  100  Oil 

Oil  000  Oil 

100  111  Oil 

101  Oil  Oil 

110  101  Oil 

111  001  Oil 
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